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CONSOLIDATION OP DATA BASE FOR ARMY GENERALIZED MISSILE MODEL 


SUMMARY 

A computer program has been written to consolidate the 
results from a series of tests on the Army Generalized Missile 
model and make them readily accessible for plotting, listing, 
and evaluation. The program is written in FORTRAN and will 
run on an ordinary minicomputer. It has the capability of 
retrieving any coefficient from the existing DATAMAN tapes and 
displaying it in tabular or plotted form. Comparisons of data 
taken in several wind tunnels and of data with the predictions 
of Program M1SSILE2 are also presented. 


INTRODUCTION 

Over the past several years a generalized missile model, 
designed and built for the Army Missile Command (MICOM), has 
been tested in several wind tunnels under a wide range of 
conditions. Some of the results of these tests have been 
recorded using Chrysler Corporation's DATAMAN system which 
produces formatted magnetic tapes. 

The purpose of the work reported here was to develop a 
computer program which would consolidate the results from these 
tests and several others (references 1-8) and make them readily 
accessible for plotting, listing, and evaluation. The program 
is written in FORTRAN and will run on an ordinary minicomputer. 
It has the capability of retrieving any coefficient from the 
existing DATAMAN tapes and displaying it in tabular or plotted 
form. 


The next eection presents the considerations involved in 
designing the data consolidation program. Following that is a 
description of the program structure and its method of 
operation. Next a user's guide to the program is given which is 
followed by a set of comparisons of data taken in different wind 
tunnels and comparisons of data with the missile aerodynamics 
prediction program, MISSILE2 (reference 9). The report is 
concluded with an appendix which contains a listing of the 
program, a description of each of the subroutines, a dictionary 
of common block variables and the program flow charts. 

PROGRAM DESIGN CONSIDERATION 

The purpose of this project is to create a computer program 
which will read a certain set of magnetic tapes and, according 
to instructions from the user, produce plots and tabulated 
listings of data from those tapes. The magnetic tapes referred 
to here are the Chrysler DATAMAN tapes listed in Table 1. The 
program is adaptable to a typical minicomputer system and uses 
existing plotting routines and hardware. 

One of the major factors in the design of the program, 
DATAMANIA, is the very large amount of data stored on magnetic 
tape. The main problem encountered in retrieving information 
from the DATAMAN tapes is that magnetic tapes are sequential 
access devices. This means that in order to read the second 
record on a magnetic tape, one must first read the first record. 
And, in order to read the last record on a tape, all preceding 
records must be read. If all the data on the tape is to be 
processed, this restriction is unimportant. However, if the 
user is searching for a particular piece of information, such as 
a dataset name, the problem becomes magnified. As each record 
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is read a comparison of some sort must be made to determine 
whether the current record is the one being sought. Even at the 
speed of modern computers this process can be quite slow. To 
overcome this problem. DATAMANIA transfers selected information 
from magnetic tape to a random access disk storage device. 

A second major factor in the design of DATAMANIA is that it 
must be adaptable to a typical minicomputer system. The system 
which has been used for the present work is described below and 
those features which are clearly typical of minicomputers are 
enumerated . 

The minicomputer facility, located at NASA/Ames Research 
Center, is organized around a Digital Equipment Corporation 
(DEC) PDP-11/70 with 192K, 16-bit words of core memory. The 
PDP-11/70 runs under DEC'S multiuser, multitasking operating 
system RSX-11M. The peripheral equipment includes a Versatec 
electrostatic printer/plotter, a nine-track tape drive, and two 
1.2-million word disk drives. Access to the computer is by 
telephone using the resource code FAC11. Resident on FAC11 is 
software which allows access to the Research Center's CDC 7600, 
a large scale computer. This software and the associated 
communication interfaces which allow the PDP-11/70 to act as a 
remote job entry station for the 7600, are not typical of 
minicomputer systems. However, a typical minicomputer system 
will include a central processing unit with a limited amount of 
core storage, a tape drive, a disk storage unit, a FORTRAN 
compiler, and some type of output device capable of drawing 
plots . 

To ensure that DATAMANIA is easily adaptable to other 
computer systems the program has been constructed in a modular 
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fashion. Each subroutine performs one function and each common 
block includes one class of data. 

A final requirement for the finished program is that it 
should be easy to use. This indicates an interactive mode of 
operation in which the user controls the direction of the 
program as it is running. This allows for a certain amount of 
trial-and-error by the user along with the benefit of neazly 
immediate results. 


GENERAL STRUCTURE AND METHOD OF OPERATION 


The program DATAMANIA does one task repeatedly. The taBk is 
to produce a frame of data which can be printed and/or plotted. 
Each frame consists of one or mere sets of independent and 
dependent variable data and pertinent descriptive information. 
The user, working from a terminal, directs the production of 
each frame by responding to prompts from the program. 

To construct a frame, the program completes the following 
steps. First, a DATAMAN tape is searched for a particular data 
set. Next, a random access file is created along with an index 
to that file. Finally, the information pertinent to the current 
frame is extracted from the random access file and stored in 
plotting arrays. If more than one dataset is required to 
complete the current frame, this process is repeated. Once a 
frame is complete it can be listed at the terminal or a plot 
file can be created for later use. 
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DATAMAN Tapes 


The Chrysler DATAMAN tapes were written on a UNIVAC 1108/3G 
computer. Due to differences between equipment It Is likely 
that the tapes will not be compatible with a typical 
minicomputer. Therefore, some preprocessing of the tapes will 
be necessary. 

Each DATAMAN tape is partitioned into datasets. A DATAMAN 
dataset represents one or more wind tunnel runs. Each dataset 
includes a unique six character identifier, a description of the 
model configuration, and a list of parametric values which arc 
constant for all runs in the dataset. The data recorded in the 
wind tunnel are prer anted as a function of two independent 
variables. The tests of Burt (Refs. 1,3,4), Hemsch (Ref. 8), 
and Schwind (Ref. 7) use Mach number and angle of attack as 
independent variables while Henderson's tests (Refs. 2,5,6) use 
either Mach number and jet coefficient or angle of attack and 
jet coefficient as the two independent variables. A description 
of the FORTRAN statements used to write a DATAMAN tape are given 
in Table 2 along with a dictionary of the variables used. 

Random Access Data Base 

In order to make efficient use of the information stored in 
each dataset and to minimize the time involved in searching a 
magnetic tape, a random access disk file is created each time a 
particular dataset is used. At the same time, the name of the 
dataset is added to a list of existing disk file datasets. In 
this way it is not necessary to perform another tape search when 
the dataset is used again. 


The random access file created by DATAMANIA is a reorganised 
and abbreviated DATAMAN dataset. The term 'random access' 
implies that each record of the file can be found in about the 
same amount of time. The system accomplishes this by numbering 
each record as it is written. In order to retrieve a particular 
record it is necessary to know its record number. DATAMANIA 
makes use of this by saving the record numbers of key records in 
an index array. 

The structure of the random access file is pictured in 
figure 1 and an example of a random access dataset is shown in 
figure 2. The first segment of the file consists of one record 
which contains identifiers for the dataset, a list of 
parameters, and a series of values which describe the size of 
the dataset. In particular, the variables NMN, MAXNAA, and NDV 
appear in this record. They represent respectively, the number 
of values of first independent variable tested, the maximum 
number of second independent variable values per first 
independent variable, and the number of dependent variables 
recorded. The record number of this record serves as the first 
index to this dataset. 

The second segment of the file represents the values of 
first independent variable encountered during the test and 
consists of NMN + 1 records. The first record, whose record 
number is the second index, contains the name of the first 
independent variable and its minimum and maximum values. This 
is followed by NMN records each of which contains MAXNAA values 
of first independent variable. When the program refers to this 
dataset, all of these values are read into an array whose 
dimensions are NMN by MAXNAA. 
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Segment three, which contains second independent variable 
data, is identical in form and sise to segment two. The first 
record of this segment contains the name of the variable and its 
minimum and maximum values. The third index for this dataset 
contains the record number of this record. 

Segment four consists of only one record which contains the 
Reynolds numbers and run numbers for this dataset. 

The fifth segment contains all of the dependent variable 
data. It is made up of NDV blocks of records, each of which are 
in the same form as the independent variable data. Each block 
begins with a dependent variable name and its minimum and 
maximum values. These record numbers are used as indices to the 
dataset and are numbered five through NDV + 4. 

The dataset shown in figure 2 is taken from Burt's Canard 
Control test in the Ames 6- by 6-Foot Wind Tunnel (reference 3) 
and is named REZ048. By comparing figures 1 and 2 it may be 
seen that the two independent variables for REZ048 are MACH and 
ALPHA. The four numbers appearing at the end of the first line 
in figure 2 are the values of the variables NMN, MAXNAA, NDV, 
and NPARM. NMN equal to two means that two Mach numbers were 
tested while a value of seven for MAXNAA means that a maximum of 
sev<»n values of angle of attack were tested for each Mach 
number. NDV equal to nine implies that nine dependent variables 
were recorded for this dataset. Their names, CN, CA, CM, CY, 
CYM, CRM, CRMB, CYB, and CNC, appear in figure 2 preceding the 
corresponding arrays of dependent variable data. A value of 
nine for the variable NPARM indicates that there are nine 
parameters for this dataset. Their names and values appear in 
the second and third lines respectively of figure 2. 


7 


The indexing scheme referred to above is implemented by 
storing the appropriate record numbers in an array named INDEX. 
Each dataset can have up to fourteen indices, depending on the 
number of dependent variables. The array INDEX has dimension 
100 by 14 and can therefore serve as an index for up to 100 
datasets. The first subscript refers to a particular dataset 
and the second to a particular record of that dataset. The 
datasets are numbered sequentially as they are created. Thus, 
the number stored in INDEX (3,1) is the record number of the 
first record of the third dataset created. In order to read the 
name of the first dependent variable of that dataset, one must 
reference the record whose record number is stored in INDEX 
(3,5). A separate array, named DSLIST, is used to store the 
names of the datasets indexed by the array INDEX. 

Method of Operation 

Once a random access file has be*n created for a particular 
dataset, the program can read the file and transfer information 
pertinent to the current frame to in-core storage areas. The 
data used by the program include all of the dataset descriptors 
and three arrays of numerical data. The three arrays, named A, 
B, and C, contain first independent variable, second independent 
variable, and dependent variable data respectively. At any one 
time, this is the only information available to the program. 

Once this information is in-core it is processed according to 
user requests. To see how this is done, consider the example 
illustrated by the following figure which shows what DATAMANIA 
knows about the hypothetical dataset AAA001. For this dataset, 
the first independent variable is Mach number, the second is 
angle of attack, and the dependent variable is CN. 


Dataset name - AAA001 Parameters 

Configuration - BODY + CANARDS PHI - o . 0 

Dependent Variable - CN PT - 7.0 

NMN ■ 3 RE AT -2.1 

MAXNAA- 8 


A 


B 


C 


0.79 

0.80 

0.81 

0.80 

0.80 

0.79 

0.81 

0.80 

1.20 

1.22 

1.21 

1.21 

1.20 

1.19 

1.20 

~nr 

1.29 

1.28 

1.28 

1.30 

1.30 

1.29 

1.28 

1.28 


Mach number 


20.1 

26.1 

30.1 

36.1 

39.8 

42.7 

45.1 

49.8 

20.2 

25.8 

28.7 

35.9 

40.1 

43.1 

45.0 

50.1 

20.1 

25.5 

29.9 

36.0 

40.8 

43.4 

44.8 

50.0 


Alpha 


5.1 

6.3 

7.4 

8.6 

9.3 

10.1 

11.3 

12.6 

6.2 

7.5 

8.1 

9.2 

10.2 

11.4 

12.7 

13.1 

7.9 

8.3 

9.4 

10.5 

11.7 

12.9 

13.6 

14.8 


If a plot of CN versus ALPHA is requested, the program will 
ask the user for a Mach number, or series of Mach numbers to be 
plotted. Suppose the user responds with Mach number ■ 0.8. 
DATAMANIA proceeds by searching the first column of A for a 
number equal to 0.8 + 0.04. The first element satisfies this 
criterion so the program transfers the first row of B and the 
first row of C to plotting arrays X and Y respectively. Notice 
that DATAMANIA has assumed that every element of the first row 
of A is near 0.8. 

To produce a plot of CN as a function of Mach number, the 
user must input a value of angle of attack as a parameter. 
Suppose a value of 30.0 is input. In this case the goal of the 
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program is to put Mach number values in array X and 
corresponding values of CN in array Y. To achieve th s end, the 
program begins by scanning the first row of B for Values of 
angle of attack that bracket 30.0; in this case, the second and 
third elements of row one. Using these numbers and the 
corresponding elements from the array C, a linear interpolation 
is performed to find the value of CN for an angle of attack of 
30.0. This value becomes the first element of the Y array. The 
corresponding Mach number is taken to be the average between the 
second and third numbers of row one in the Mach number array. 
Again, it is assumed that the Mach numbers will not vary 
significantly. DATAMANIA then repeats this process for each 
remaining row. 

Since each dataset represents a test performed at one roll 
angle, PHI, it is necessary to use a series of datasets to 
produce a plot with roll angle as the independent variable. 
DATAMANIA depends on the user's good judgment in selecting these 
datasets. To produce a meaningful plot, the parameters from 
each dataset should be identical (excepting PHI, of course), as 
should the configurations. DATAMANIA processes these datasets 
one at a time, performing the same task each time. First, the 
parameter list is searched for the word PHI. The corresponding 
parametric value is placed in the X array. To find the value cf 
CN which is associated with this roll angle, bilinear 
interpolation is used. The user will have specified a Mach 
number and an angle of attack for which a CN value is desired. 
The Mach number array and the angle of attack array are searched 
for bracketing values from which CN is to be determined. Then, 
CN is determined by first interpolating in Mach number and then 
in angle of attack. The resulting CN is stored in the Y array 
for later plotting or listing. 
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Plotting 


An important aspect of the design of DATAMAN1A is its 
adaptability to a variety of plotting packages. This was 
accomplished by passing all plot information in two common 
blocks; one containing the arrays to be plotted and the other 
containing textual data which is to appear on the plot. 

For the present implementation of DATAMANIA the PLOTA/PLOTI 
(reference 10) routines were chosen. This package, developed at 
Nielsen Engineering & Research, Inc. (NEAR), was selected for 
several reasons. First, it is simple to use; a single FORTRAN 
CALL statement will generate an X,Y plot. Second, automatic 
scaling or user defined scales are available. Third, the code 
is readily available and the necessary host routines are aJ ready 
resident on the computer system which was used. Finally, the 
entire package is clearly and thoroughly documented. 

Output from the PLOTA/PLOTI routines consists of plot 
information in the Simple Graphics Package (SGP) File format 
(reference 11). This intermediate file, which can be read and 
edited, can then be disposed to a device-specific post processor 
to create CRT displays or hard copy plots. Post processors have 
been written (reference 11) for the following devices: 

1 ) Evans and Sutherland Picture System 

2) Versatec Electrostatic Printer/Plotter 

3) IMLAC - PDS1 Graphic Display Scope 

4) Zeta Pen Plotter 

5 ) SC4020 Hard Copy Device 

The plots appearing in this report were originally produced 
on a Versatec printer/plotter. 
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Implementation on DEC PDP-11/70 


The implementation of DATAMANIA on the DEC PDP-11/70 
described above is governed by several factors. First, NEAR is 
physically separated from the computer by more than a mile. The 
system does not have a full-time operator to respond to requests 
for tape mounts or other auxiliary services. Also, the system 
supports many users with the result that public disk storage 
space is limited. Finally, the DATAMAN tapes created by 
Chrysler Corporation cannot be read by DEC machinery. Because 
of these limitations, it was decided to make use of the 
communicating link between the minicomputer and the large scale 
CDC 7600. 

The first necessary step for running DATAMANIA on the 
PDP-11/70 is to translate the DATAMAN data to a format 
compatible with the DEC system. This is accomplished by reading 
the tapes on the 7600 according to the FORMAT statements in 
Table 2 and then writing the data to a 7600 disk file. In this 
process the arrays NAMDM, DM, and GO, which are not used on any 
of the DATAMAN tapes, are deleted from each dataset. The disk 
file created in this way can be sent to the PDP-11/70 using the 
communicating link between the two computers. Once this file is 
captured on the DEC system it is given the name FOR001.DAT. By 
DEC convention, this name is associated with logical unit l for 
FORTRAN programs; i.e., when DATAMANIA reads data from logical 
unit 1, it reads the file named FOR001.DAT. 

The fundamental executable unit on the DEC minicomputer is a 
task. To make DATAMANIA a task, the source code, as shown in 
the appendix, is entered in the system as a text file by using 
one of the system's editors. Next, the text file is translated 
into an object module by using the FORTRAN compiler. The object 
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module is then converted to a task image using DEC'S Task 
Builder program. The program DATAMANIA is then ready to run. 

The only nonstandard FORTRAN language features used in DATA- 
MANIA concern the creation and use of a random access file. A 
random access file is created using the statement OPEN with 
appropriate parameters. Reading the file requires the use of a 
READ statement whose syntax is nonstandard. Implementation of 
this program on another type of machine will involve replacing 
each occurrence of these two statements with their machine 
dependent counterparts. 

The PLOTA/PLOTI plot pack*».?e used by DATAMANIA writes its 
output to logical unit sever.. On the PDP-ll/70 this means that 
a disk file named FOR007.DAT is created. A standard plot 
created by the PLOTA/PLOTI package has a bottom axis and a left 
side axis only. The plotting routines themselves have been left 
in standard form except for the addition of two subroutines, 

RITAX and TOPAX, which draw the top axis and right side axis. 

At the end of a DATAMANIA session, plots are generated by 
disposing file FOR007.DAT to the appropriate SGP post-processor. 

USER'S GUIDE 

Success in the use of DATAMANIA depends to a large degree on 
the preparation of the user. Since the program is intended to be 
a model for similar implementations on other computer systems, it 
was considered important to keep the code simple. As a result, 
many error detecting and processing devices have been omitted 
which would be valuable in a version oriented towards production. 

To prepare for a DATAMANIA session the user should make a 
list of the plots to be generated and a list of the datsets 
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necessary to make those plots. This Information can be compiled 
from the wind tunnel reports, references 2 through 8. Once this 
is done, and the necessary DATAMAN tapes are ready for mounting, 
the user is ready to use DATAMANIA. 

DATAMANIA builds frames of data by asking the user for 
descriptive information. The first prompt to the user is for 
independent variable name. He responds by typing either ALPHA, 
MACH, CRT, or PHI. This response determines the program's mode 
of operation until the current frame is complete. 

If either ALPHA, MACH, or CRT are named as the independent 
variable, DATAMANIA will ask for the name of the dataset to be 
used and the name of the dependent variable being plotted. 

Next, the user is asked to input a series of parametric values. 
The parameter used will depend on the independent variable and 
the dataset being used. One curve will be produced for each 
parameter. When this data has been assembled the user has the 
option of putting more curves on the same frame using either the 
same dataset or another dataset. A maximum of five curves per 
frame is allowed. 

When a frame is to be generated using PHI as the independent 
variable, the user must input a list of dataset names to be used. 
It is the responsibility of the user to ensure that these data- 
sets are compatible with one another and that each represents a 
test performed at a different roll angle. After the dataset list 
is input, the program prompts for a dependent variable name. Now, 
the user has the option of producing either a series of curves 
each representing a different Mach number and one angle of attack 
or representing different alphas and one Mach number. DATAMANIA 
determines this by printing "ENTER NAME OF PARAMETER TO BE HELD 
CONSTANT (MACH OR ALPHA) - ". 
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As each frame is completed DATAMANIA asks the user if a 
tabulation is desired. If the response is "YES", the contents 
of the current frame are listed at the terminal. Next, the user 
is asked if a plot file should be created. If so, he is 
prompted for the type of scaling to be used. The scaling 
options are default scaling, same-as- last-plot scaling, and user 
defined scaling. The default scaling option will produce plots 
which include all data points within the smallest range which 
will allow reasonable labeling of the tick marks. To generate 
plots with user defined scales the user must input minimum and 
maximum values for both the independent variable and the 
dependent variable. 

Sample Session 

The following sample session with DATAMANIA will illustrate 
its use in generating three frames of data each using a 
different independent variable. Underlining has been used to 
show user responses. All three frames are tabulated as well as 
plotted . 


> RUN DATAMANIA 


********** WELCOME TO DATAMANIA ********** 


DO YOU WANT INSTRUCTIONS? YES 


DATAMANIA IS DESIGNED TO READ DATAMAN TAPES AND 
PRODUCE PLOTS OF COEFFICIENT DATA AS A FUNCTION OF 
EITHER ANGLE OF ATTACK (ALPHA), MACH NUMBER, THRUST 
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COEFFICIENT (CRT), OR ROLL ANGLE (PHI). IN ORDER 
TO DO THIS, YOU NEED TO KNOW THE NAMES OF THE 
DATASETS WHICH CONTAIN THE INFORMATION TO BE PLOTTED. 

AT THE BEGINNING OF EACH FRAME, DATAMANIA WILL 
PROMPT FOR THE INDEPENDENT VARIABLE NAME. YOUR 
RESPONSE OF ALPHA, MACH, CRT, OR PHI WILL DETERMINE 
THE MODE OF PROGRAM OPERATION. FOLLOWING THIS YOU 
WILL BE PROMPTED FOR INFORMATION NEEDED TO COMPLETE 
THE FRAME. 

AS EACH FRAME IS COMPLETED YOU WILL BE ASKED 
FIRST IF YOU WANT A TABULATION OF THE CURRENT 
FRAME, AND THEN, IF YOU WANT TO CREATE A PLOT FILE 
OF THE FRAME. 


ENTER INDEPENDENT VARIABLE NAME - ALPHA 


ENTER DATASET IDENTIFIER - RAW002 


DO YOU WANT A SUMMARY OF RAW002? YES 
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SUMMARY OF RAW002 


CREATED 22 SEP 77 
CONFIGURATION - BODY + TAILS 

PARAMETERS 


RE/FT * 2.100 
PHI * 0 . 000 
PT * 7 . 000 

MACH NUMBERS 


INDEX RANGE 

1 FROM 0.790 TO 0.790 

2 FROM 1.210 TO 1.220 

3 FROM 1.280 TO 1.290 

DEPENDENT VARIABLES 


CN 

CY 

CA 

CLM 

CYN 

CBL 

PHI-M 

PHI-C 

PHI-T 

BETA 


ENTER DEPENDENT VARIABLE NAME - CN 

ENTER NUMBER OF MACH NUMBERS TO BE PLOTTED - 2 

ENTER 2 MACH NUMBERS (SEPARATE WITH COMMAS) - .8,1.2 

ANY MORE DATA ON THIS FRAME? NO 

TABULATE DATA? YES 


$$$$$$$$$$$$$ TABULATION OF CURRENT FRAME 

$ 

$ 

? 

$ DATASET - RAW002 
$ CONFIGURATION - BODY + TAILS 
$ PARAMETRIC VALUES 

$ 

$ RE/FT 2.100 

$ PHI 0.000 

$ PT 7 . 000 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 


18 


FOR MACH >0.8 FROM RAW002 


$ 


$ 

ALPHA 

CN 

$ 

20.0700 

4.3457 

$ 

21.9700 

4.7439 

$ 

23.8400 

5.2159 

$ 

26.8300 

6.0867 

$ 

29.7400 

6.8728 

$ 

32.7600 

7.6546 

$ 

34.6600 

8.1298 

$ 

38.5600 

9.2253 

$ 

41.6400 

10.0630 

$ 

44.7200 

11.0090 

$ 

47.5600 

11.9860 

$ 

49.2400 

12.4710 


$ 

$ 

$ ALPHA CN FOR MACH -1.2 FROM RAW002 

$ 20.1300 5.2252 

$ 22.0600 5.8659 

$ 23.9500 6.4432 

$ 26.9800 7.5335 

$ 29.9400 8.6510 

$ 33.0000 9.9209 

$ 34.9300 10.7130 

$ 38.8800 12.2730 

$ 41.9800 13.3800 

$ 45.0600 14.3030 

$ 47.8900 14.9450 

$ 49.5700 15.2940 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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■ 

> 

* 


CREATE PLOT FILE? YES 

USE DEFAULT SCALING? YES 

WANT TO MAKE ANOTHER FRAME? YES 

ENTER INDEPENDENT VARIABLE NAME - MACH 

ENTER DATASET IDENTIFIER - RAW002 

DO YOU WANT A SUMMARY OF RAW002? NO 

ENTER DEPENDENT VARIABLE NAME - CN 

ENTER 3 ALPHAS (SEPARATE WITH COMMAS) - 20. ,30, ,50. 
ANY MORE DATA ON THIS FRAME? NO 
TABULATE DATA? YES 


$$$$$$$$$$$$$ TABULATION OF CURRENT FRAME 

$ 

§ 

$ 

$ DATA SET - RAW002 
$ CONFIGURATION - BODY + TAILS 
$ PARAMETRIC VALUES 

§ 

$ RE/FT 2.100 

$ PHI 0.000 

$ PT 7 . 000 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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$ 

$ 


$ 

MACH 

CN 

FOR 

ALPHA-20.0 

FROM 

n 

o 

0 

1 

$ 

0.7900 

4.3310 





$ 

1.2150 

5.1820 





$ 

1.2900 

5.1057 





$ 







$ 







$ 

MACH 

CN 

FOR 

ALPHA-30.0 

FROM 

RAW002 

$ 

0.7900 

6.9 





$ 

1.2150 

8.6759 





$ 

1.2900 

8.6621 






$ 

$ 


§ 

MACH 

CN 

FOR ALPHA-50.0 FROM RAW002 

$ 

0.7900 

12.6904 


$ 

1.2100 

15.3833 


$ 

1.2850 

15.5535 



$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

CREATE PLOT FILE? YES 

USE DEFAULT SCALING? YES 

WANT TO MAKE ANOTHER FRAME? YES 

ENTER INDEPENDENT VARIABLE NAME - PHI 

ENTER NUMBER OF PHI VALUES TO BE PLOTTED - 4 
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ENTER 4 DATASET NAMES 

EACH FOLLOWED BY CARRIAGE RETURN 

RAW002 

RAW003 

RAW004 

RAW005 

ENTER DEPENDENT VARIABLE NAME - CN 

ENTER PARAMETER TO BE HELD CONSTANT (MACH OR ALPHA) - MACH 
ENTER VALUE OF CONSTANT PARAMETER - jJ) 

ENTER NUMBER OF ALPHA VALUES TO BE PLOTTED - 3 
ENTER 3 ALPHA'S (SEPARATE WITH COMMAS) - 20. ,30. ,50, 

DO YOU WANT A SUMMARY OF RAW002? NO 
DO YOU WANT A SUMMARY OF RAW003? NO 
DO YOU WANT A SUMMARY OF RAW004? NO 
DO YOU WANT A SUMMARY OF RAW005? NO 
TABULATE DATA? YES 
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$$$$$$$$$$$$ TABULATION OF CURRENT FRAME 

$ 

$ USINO DATASETS 
$ RAW002 

$ RAW003 

$ RAW004 

$ RAW005 

$ 

$ CONFIGURATION - BODY ♦ TAILS 

$ 

$ PARAMETRIC VALUES 

$ 

$ RE/FT 2.100 

$ PT 0.000 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 


MACH - 0.8 


$ 

$ 

$ PHI 

$ 0.0000 

$ 10.0000 

$ 30.0000 

$ 45.0000 

$ 

$ 

$ PHI 

$ 0 . 0000 

$ 10.0000 

$ 30.0000 

$ 45.0000 

$ 


CN FOR ALPHA-20.0 

4.3513 
4.3836 
4.1957 
4.2816 


CN FOR ALPHA-30.0 

6.9814 
6.8639 
6.7752 
6.5203 


23 


FOR ALPHA-50.0 


$ 

$ PHI CN 

$ 0.0000 12.7545 

$ 10.0000 -5.6740 

$ 30.0000 11.8958 

$ 45.0000 11.5471 

$ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

CREATE PLOT FILE? YES 


USE DEFAULT SCALING? YES 
WANT TO MAKE ANOTHER FRAME? NO 
********** END OF SESSION **** 


DATA COMPARISONS 

The consolidated data base for the Army Generalized Missile 
can be divided into three categories: 

1. Henderson's plume interaction tests, references 2,5,6. 

2. Burt's tests of nose mounted canard configurations, 
references 1,3,4. 

3. NEAR high ancile-of-attack tests, references 7,8. 

Comparisons of data taken in different wind tunnels can be made 
for the first and third categories and are presented below. 
Also, comparisons of data and the predictions of Program 
MISSILE2 (reference 9) are presented below for some of the 
tests in categories 2 and 3. 


& : 


24 


Comparisons of Data Taken in Different Wind Tunnels for 
Plume Interaction Tests 

A comprehensive series of plume interaction tests were run in 
the Calspan 8-Foot Transonic Wind Tunnel and the AEDC 16T and 
16S Wind Tunnels for several body-alone and body-tail configura- 
tions. Seme of the test parameters at the two facilities were 
identical and provide a means of checking the sensitivity of the 
results to test conditions. Figure 3 provides geometrical 
details of the Army Generalized Missile model. Comparisons are 
presented for the following configurations. 

(a) Body-alone with Army nose for the following conditions: 


(b) 


Mach 

Radial Thrust 
Coefficient/ CRT 

Figure 

0.4 

0/ 12 

4a 

1.0 

0, 6 

4b 

1.25 

6 

4c 


Body with Army 

nose and Army FI tail 

fins mounted in the 

forward hinge- 

line position for the 

following conditions 

Mach 

Radial Thrust 
Coefficient/ CRT 

Figure 

0.4 

0 

5a 

1.0 

0, 12 

5b 

1.25 

6 

5c 
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(c) Body with Army nose and Army F2 tail fins mounted in the 
forward hinge-lit.e position for the following conditions i 


Radial Thrust 

Mach Coefficient# CRT Figure 


1.0 


0 # 6 


6 


The pitching-moment coefficient is shown for comparison 
since it is usually a more sensitive variable than normal force. 
In general, except for zero shifts at M*, ■ 0.4 for the body 
alone and * 1.25 for the Fin 1 case, the comparisons for all 
three configurations are good. 


Comparisons of Data From Nose-Mounted Canard-Configuration 
Tests with Predictions of Program MISSILE2 

In this section, comparisons of data with the prediction of 
Program MISSILE2 (reference 9) are presented for a particular 
canard configuration with inline and with interdigitated tail 
fins. The comparisons are for three Mach numbers (M^ ■ 0.9, 
1.25, 2.00) and for roll angles of 0 and 45 degrees. The 
configuration consists of the Army Generalized Missile body with 
the Army nose, the Army canards mounted at the aft canard 
hinge-line position, and the Army FI tail fins mounted at the 
aft tail hinge-line position. 

Results for overall normal-force coefficient with zero roll 
angle and no fin deflection are presented in figure 7. While 
the Mach number and angle of attack trends are correct, the 
theory underpredicts the data by about 10-15%. The rest of the 
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comparisons for this configuration are for the pitching-moment 
coefficient which is a sensitive indicator of the ability of the 
theory to predict longitudinal stability. The results are 
presented in figures 8, 9 and 10 for ■ 0.9, 1.25 and 2.0 
respectively. (Note the difference in moment center for the M* 

■ 2.0 case.) Except for the 15° canard fin deflection case, 
the agreement is good for ■ 0.9 and 1.25. The worst 
agreement (excluding the 15° fin deflection cases) gives an 
error of roughly 0.4 body diameters in the position of the 
longitudinal center of pressure. For ■ 2.0# the maximum 
error is nearly one body diameter. 

Comparisons of Data from NEAR High Angle-of-Attack Tests 
With Prediction of Program MISS1LE2 

In this section, comparisons of data with the predictions of 
Program MISSILE2 (reference 9) are presented for another canard 
configuration with inline tail fins. The Navy nose, canards and 
tail fins shown in figure 3 were used. The comparisons are for 
subsonic and supersonic speeds (M^ * 0.8, 1.3), several roll 
angles and with and without fin deflection. Since Program 
MISSILE is especially useful for the determination of lateral 
stability and control-cross coupling effects, comparisons are 
presented for the induced rolling-moment coefficient. 

Results for zero fin deflection at roll angles of 10° 
and 300 are presented in figure 11. Note that the data 
taken in the NASA/ARC 6- by 6-Foot Wind Tunnel and NASA/ARC 
11 -Foot Wind Tunnel overlap for angles of attack between 20 
and 24°. The agreement for the two tunnels is satisfactory. 

The agreement for the roll angle case of 30° is excellent 
for angles of attack up to 40° even to the point of repro- 
ducing the Mach number effects between angles of attack of 
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150 an d 350. For IOO roll, the agreement is acceptable 
for angles of attack up to 40 ® for purposes of design. 

Results for 150 deflection of the yaw control fins at roll 
angles of Oo and 300 are presented in figure 12 . Again, agree- 
ment of data from the two wind tunnels is satisfactory. For 
this case, the theory predicts the overall level and trends of 
induced roll adequately for design purposes. However, definite 
improvement of the program is needed for the higher angle of 
attack range. 
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APPENDIX 


Description of Subroutines 

MAIN - This is the driver for the rest of the program. It 

begins by opening the random access file (logical unit 3) 
and initializing counters and flags. Then, for each frame, 
MAIN asks the user for the independent variable name and 
the variable MODE is set. Next, subroutine MYSTRO is 
called to produce a frame. When the frame is complete, 
control returns to MAIN which then responds to user 
instructions to tabulate and/or plot the frame. See flow 
chart in figure 13. 

MYSTRO - This subroutine does most of the interaction with the 
user. Its job is to obtain a description of the desired 
frame and to direct its production. See flow chart in 
figure 14. 

CREATO - Given a dataset name, this subroutine first searches a 
list of existing datasets for that name. If the dataset 
name appears on the list, the variable DSNUM is given an 
appropriate value and control is returned to MYSTRO. If 
the random access dataset does not exist, a DATAMAN tape is 
searched, a new random access dataset is created, and its 
name is added to the list of existing dataset names. See 
flow chart in figure 15. 

FINDER - This subroutine reads the random access file created 
by CREATO. Given a dependent variable name, FINDER fills 
the arrays A, B, and C with data from the current dataset. 
Array A contains first independent variable data, B 
contains second independent variable data, and C is an 
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array of dependent variable data. This routine also 
retrieves identifying information for the dataset such as 
the configuration , date of creation, and the names and 
values of the parameters. See flow chart in figure 16. 

ABC - Subroutine ABC processes the arrays A, B, and C which are 
provided by subroutine FINDER. Depending on the 
independent variable of the current frame, this routine 
writes information into arrays X and Y which repzesent the 
abscissa and ordinate values to be plotted. Linear 
interpolation is done between the discrete elements of A, 

B, and C. See flow chart in figure 17. 

SORTOR - This subroutine is called only if roll angle is the 
independent variable. Its purpose is to put the plotting 
ordinates in ascending order. 

TABULA - This routine directs the formatting and listing of the 
current frame at the terminal. 

PLOTOR - Subroutine PLOTOR is called if a plot file is to be 
created. After asking the user for scaling information, 
it directs the production cf an SGP plot file. 


Dictionary of Common Block Variables 

Common block BETA contains information concerning the 
dataset currently being processed. Its variables are defined in 
subroutine FINDER and processed by subroutine ABC. 

A an array of first independent variables 

B an array of second independent variables 
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c 

NBP 

SETNAM 

CONFIG 

DVNAM 

NMN 

MAXNAA 

NPARM 

NAMPAR 

PVAL 

FEEVAL 

TMACH 

TALPH 

Common 

MODE 


an array of data for a specific dependent variable 

an array contaning the number of entries recorded 
in each row of A, B, and C 

name of current dataset 

configuration of current dataset 

name of dependent variable under consideration 

number of runs in dataset 

largest value in array NBP 

number of parameters in current dataset 

names of parameters 

parametric values 

value of roll angle in parameter list 

Mach number under consideration 

angle of attack under consideration 

block COUNT contains counters and flags. 

flag indicating independent variable for current 
frame (1 ■ first independent variable, 2 * second 
independent variable, and 3 ■ roll angle) 


31 


NCURV 

number 

of 

NSETS 

number 

of 

NAPARM 

number 

of 

NFEE 

number 

of 


MODP 


NPT 


FLAG 


number of datasets used to make current frame 


independent variable 

flag indicating active parameter when roll angle is 
the independent variable (1 ■ first independent 
variable, 2 ■ second independent variable) 

number of points plotted when roll angle is the 
independent variable 

flag indicating an error condition 


Common block CURVS contains information defining the curves 
to be plotted. Its variables are set in subroutine ABC and used 
by PLOTOR. 

X array of abscissae to be plotted 

Y array of ordinates to be plotted 

NP number of points in each curve 

AVAL values of the active parameter for each curve 

CPAR value of parameter held constant 
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Common block FILE3 appears in all subroutines that read or 
write to the random access data base, logical unit 3. 

NREC3 the variable associated with unit 3. It is 

incremented by one after each read or write 
operation on unit 3 

NEXREC after each data base is written on unit 3, this 

variable is set to the next available record number 

Common block IDENTS contains identifiers for the frame 
currently being constructed. 

SETID the name of dataset input by the user 

VAR contains the names of both independent variables 

for the current dataset 

IVID independent variable name for the current frame 

DVID dependent variable name input by the user 

ACPARM values of the active parameters input by the user 

NAMES the names of the datasets to be used when roll 

angle is the independent variable 

Common block NDEX contains data used to index the random 
access datasets which are created by subroutine CREATO. 

INDEX an array of indices for the random access data 

base 
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DSLIST a list of existing random access datasets 

DSNUM the identifying number of the dataset currently 

being used 

SETCNT the number of datasets in D8LI8T 

Common block TEXT contains textual identifiers which appear 

on plots and tabulations of the current frame. 

IDCURV an array containing the dependent variable name 

and dataset name associated with each curve in the 
current frame 

IDSETS a list of datasets used in the current frame 

IDCFIG a list of configurations associated with the data- 

sets in IDSETS 

IDPARMS an array containing a list of parameters for each 

dataset in IDSETS 

NPARMS the number of entire in IDPARMS 

PARVAL the parametric values for each parameter in IDPARMS 

Program Listing 
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1 


PROGRAM DATAMANIA 
540/C 

IMPLICIT INTEGER (A-Z) 

COMMON /COUNT/ MODE , NCURV , N8ETS , NAPARM , NT EE , MODP , NPT , FLAG , 

* IVFLAG, ACFLAG, MACP 

COMMON /IDENTS/ SETID(3 ) , VAR( 2 , 3 ) , IVID( 3 ) , DVID( 3 ) , ACPARM( 5 ) , 

* NAMES (iO, 3) 

COMMON /FILE3/ NREC3 , NEXREC 

COMMON /NDEX/ INDEX (14, 100) , DSLIST(100, 3) , D8NUM, SETCNT 

DATA YES /1HY/ , 

* MOCK /2HMA/ , 

* ALFA / 2HAL/ , 

* FEE /2HPH/, 

* CRT /2HCR/ , 

* ZERO /0/ 

OPEN (UNIT-3 , TYPE- ' NEW' , ACCESS- ' DIRECT ' , FORM- ' UNFORMATTED ' , 

* RECORDS IZE-64, INITIALSIZE-100, DISPOSE-' DELETE' , 

* ASSOC IATEVARIABLE-NREC3 ) 

NEXREC— 1 
SETCNT— 0 
FLAG-0 
DSNUM-0 

WRITE DATAMANIA HEADER AND PROMPT FOR INSTRUCTIONS 

WRITE (5, 5001) 

READ (5,5002) ANSWER 

IF ( ANSWER. EQ. YES) CALL NSTRUC 

PROMPT FOR INDEPENDENT VARIABLE 

50 NSETS-0 
NCURV— 0 
MODE-O 

WRITE (5, 5003) 

READ (5,5009) IVID 

IF (IVID(l) .EQ.ALFA .OR. IVID(1 ) . EQ.MOCK .OR. 

* IVID(1 ) . EQ.CRT .OR. IVID(1 ) .EQ.FEE) GO TO 100 


UNACCEPTABLE IND. VAR. - WRITE MESSAGE, TRY AGAIN 

WRITE (5, 5004) 

GO TO 50 

100 CALL MYSTRO 

IF ( FLAG. EQ. ZERO) GO TO 150 
FLAG-0 

WRITE(5, 5010) 

GO TO 50 


I. 



PROMPT FOR TABULATION 
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C 

150 CONTINUE 

WRITE (5, 5005) 

READ (5,5002) ANSWER 
IF ( ANSWER. EQ. YES) CALL TABULA 

PROMPT FOR PLOT FILE 

WRITE (5, 5006) 

READ (5,5002) ANSWER 
IF ( ANSWER. EQ. YES) CALL PLOTOR 

PROMPT FOR MORE 

WRITE (5, 5007) 

READ (5,5002) ANSWER 
IF (ANSWER. EQ. YES) 00 TO 50 

END OF 8ESSION 

WRITE (5, 5008) 

STOP 


FORMATS******************* ******* ****** 

5001 FORMAT ( / ///IX, 10(1H*) ,22H WELCOME TO DATAMANIA ,10(1H*)//// 

* 27H$DO YOU WANT INSTRUCTIONS? ) 

5002 FORMAT (Al) 

5003 FORMAT (/35H$ENTER INDEPENDENT VARIABLE NAME - ) 

5004 FORMAT (/43H INDEPENDENT VARIABLE NAME IS UNACCEPTABLE/ 

*59H THE FOUR ACCEPTABLE IND. VARIABLES ARE ALPHA, MACH NUMBER,/ 
*59H CRT, AND PHI. DATAMANIA READS ONLY THE FIRST TWO LETTERS/ 
*54H OF YOUR RESPONSE SO YOU MAY ENTER AL, MA, CR, OR PH.) 

5005 FORMAT ( / 1 6H ^TABULATE DATA? ) 

5006 FORMAT (/19H$CREATE PLOT FILE? ) 

5007 FORMAT ( /29H$WANT TO MAKE ANOTHER FRAME? ) 

5008 FORMAT (//IX, 13 (1H*),16H END OF SESSION ,13(1H*)//) 

5009 FORMAT (3A2) 

5010 FORMAT (//45H IND. VAR. NAME UNACCEPTABLE FOR THIS DATASET/) 


END 


nnn n n n n n ooooo ooooo 


SUBROUTINE ABC 


THIS SUBROUTINE EXTRACTS DATA FROM THE ARRAY8 
GENERATED BY FINDER AND PREPARES IT FOR USE 
BY THE PLOTTING AND TABULATING ROUTINES 

COMMON /COUNT/ MODE , NCURV , NSETS , NAPARM , NFEE , MODP , NPT , FLAG , 

* 1VFLAG , ACFLAG , MACP 
COMMON /BETA/ A(10, 20) ,B(10, 20) ,C(10, 20) ,NBP(10) , 

* SETNAM( 3 ) , CONFIG ( 24) , DVNAM( 3 ) , NMN, MAXNAA, NPARM, 

* NAMPAR (3,10), PVAL (10), FEEVAL, TMACH , TALPH 
COMMON /IDENTS/ SETID(3 ) , VAR(2 , 3 ) , IVID( 3 ) ,DVID(3 ) , ACPARMt 5 ) , 

* NAMES (10, 3) 

COMMON /CURVS/ X ( 20, 5 ) , Y( 20, 5 ) , NP( 5 ) , AVAL(5 ) ,CPAR 
COMMON /TEXT/ IDCURV(6, 5) , IDSETS(3 , 10 ) , IDCFIG(24, 5 ) , 

* IDPARM(3, 10, 5) ,NPARMS(5) ,PARVAL(10, 5) 

INTEGER SETNAM, CONFIG, DVNAM, SETID, DVID, ACFLAG, VAR 
YLINE (OR1 , AB1 , OR2 , AB2 , ORP )-ABl+ ( AB2-AB1 ) * (ORP-OR1 ) / (OR2-OR1 ) 
DATA EPS /0.04/ 

NPARMS ( NSETS ) -NPARM 
DO 10 1-1,3 

I DSETS ( I , NSETS ) -SETNAM ( I ) 

10 CONTINUE 

DO 20 1-1,24 

IDCFIG ( I , NSETS ) -CONFIG ( I ) 

20 CONTINUE 

DO 30 M-l. NPARM 
PARVAL ( M , NSETS ) -PVAL ( M ) 

DO 30 1-1,3 

IDPARM( I , M, NSETS ) -NAMPAR ( I , M) 

30 CONTINUE 

IF ( MODE . EQ . 3 ) GO TO 500 

DO 400 N-l , NAPARM 
NCURV— NCURV+1 
AVAL (NCURV ) — ACPARM ( N ) 

DO 40 1-1,3 

IDCURV ( I , NCURV ) -DVNAM ( I ) 

IDCURV( 1+3 , NCURV) -SETNAM ( I ) 

40 CONTINUE 

USE 2ND IND. VAR. 

IF ( MODE . EQ . 1 ) GO TO 200 
C 

DO 50 I-l.NMN 

IF (ABS(ACPARM(N)-A(I, 1) ) .LT.EPS) GO TO 100 
50 CONTINUE 
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noon n 00000 


1 - 1-1 

100 INC1-I 

NP (NCURV ) *NBP ( I NCI ) 

DO 150 J*1 , NBP ( INC1 ) 

X ( J , NCURV ) *B ( INC1 1 J ) 

Y(J,NCURV)-C{INC1, J) 

150 CONTINUE 

GO TO 400 

USE 1ST IND. VAR. 

200 CONTINUE 

NP( NCURV )=NMN 

DO 330 M*1 , NMN 

DO 300 1*2 , MAXNAA 

II=I 

IF ( ACPARM ( N ) . LE . B ( M , I ) ) GO TO 350 
300 CONTINUE 

350 1=11-1 

X (M, NCURV)* (A(M, II )+A(M, I))*0.5 

Y(M,NCURV)-YLINE(B(M,I),C(M,I) ,B(M, II ) ,C(M, II ) , ACPARM(N) ) 
330 CONTINUE 
400 CONTINUE 

RETURN 

IND. VAR. IS PHI 

500 CONTINUE 

DO 750 N=1 , NAPARM 

NCURV*N 

AVAL ( NCURV ) =ACPARM ( N ) 

DO 530 1*1,3 

IDCURV ( I , NCURV ) =DVNAM ( I ) 

530 CONTINUE 

IF (MODP.EQ.l) TALPH=ACPARM(N) 

IF (MODP.EQ.2) TMACH*ACPARM(N) 

DO 550 1*2, NMN 

II*I 

IF ( TMACH . LE . A ( I , 1 ) ) GO TO 600 
550 CONTINUE 
C 

600 1=11-1 
C 

DO 650 J*2, MAXNAA 

JJ*J 

IF (TALPH.LE.B(I, J) ) GO TO 700 
650 CONTINUE 
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c 


c 

c 

c 

c 

c 

c 


700 J-JJ-1 

X ( N PT , NCURV ) -FEEVAL 

C13-YLINE(B(I,J),C(I, J) ,B(I, JJ) ,C(I, J J ) , TALPH) 
C24*YLINE ( B( II « J ) , C( 1 1 , J ) , B( I I , JJ ) , C ( I I , JJ ) , TALPH) 

Y (NPT $ NCURV )*YLINE (A( I,1) # C13,A(II,1) #C24, TMACH) 

750 CONTINUE 

RETURN 

END 


non ooo nnn n n n o n onooo 


SUBROUTINE CREATO 


THIS ROUTINE CREATES A DATABASE AND AN INDEX 
FOR THE DATASET UNDER CONSIDERATION. 

IMPLICIT INTEGER (A-Z) 

COMMON /COUNT/ MODE, NCURV, NSETS, NAPARM, NFES. MODP, NPT, FLAG, 

* IVFLAG, ACFLAG, MACP 
COMMON /FILE3/ NREC3 , NEXREC 

COMMON /IDENTS/ SETID( 3 ) , VAR(2 , 3 ) , IVID( 3 ) , DVID( 3 ) , ACPARM( 5 ) , 

* NAMES (10,3) 

COMMON /NDEX/ INDEX (14, 100 ) , DSLIST ( 100, 3 ) , DSNUM, SETCNT 

REAL DUM( 10 ) , XLIM( 2 ) , APLIM( 2 ) , DLIM( 10, 2 ) , PARVALt 10 ) , 

* A( 10, 20 ) , ACPARM 

DIMENSION DSNAM(3) ,DATE(6) ,CONFIG(24) ,DVAR(10, 3) , 

* PARMS (10,3) , RUNNO (10,2), NAAX (10) 

DATA YES /1HY/ 

IF (SETCNT. EQ.O) GO TO 10 

SEARCH ARRAY OF EXISTING DATASETS 


DO 5 N-l , SETCNT 
DSNUM-N 

IF (SETID(1 ) .NE.DSLIST(N, 1 ) ) GO TO 5 
IF ( SETI D( 2 ) . NE . DSLIST (N , 2 ) ) GO TO 5 
IF ( SETID( 3 ) . EQ. DSLIST (N, 3 ) ) GO TO 65 
5 CONTINUE 

SEARCH FOR DATASET IDENTIFIER 

10 READ( 1,1001, END-70) DSNAM 

IF ( DSNAM (1 ) .NE.SETID(1 ) ) GO TO 10 
IF ( DSNAM ( 2 ) • NE . SETID ( 2 ) ) GO TO 10 
IF ( DSNAM ( 3 ) .NE . SETID (3 ) ) GO TO 10 

YOU’VE FOUND IT - NOW CREATE 

BACKSPACE 1 
SETCNT-SETCNT+1 
DSNUM-SETCNT 
C 

DO 20 1-1,3 

20 DSL1ST(DSNUM.I)-DSNAM(I) 

C 

READ(1, 1002) DSNAM. DATE, CONFIG. NMN, MAXNAA, NDV, NPAR 
C 

READ( 1 . 1003 ) ( (VAR ( I , J ) , J-l , 3 ) , 1-1 , 2 ) , 

* ( ( DVAR l I , J ) , J-l , 3 ) , I— 1 , NDV ) , 

* ( ( PARMS ( I , J ) , J-l , 3 ) , I— 1 , NPAR) , 

* (DUM(I), 1-1,30) 

n 

READ ( 1 , 1004 ) ( ( RUNNO ( I , J ) , 1-1 . NMN ) , J-l , 2 ) , ( NAAX ( I ) , 1-1 , NMN ) 
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c 

READ( 1 . 1005 ) (DUMl I ) , 1-1 . 2 ) , (XL1'M( I ) . 1-1 . 2 ) , ( APLIM( I ) . 1-1 , 2 ) , 

* ((DLIM(I,J),I-1,NDV),J-1,2), 

* (PARVAL(I).I-l.NPAR), (DUM( I ) , 1-1 . 7 ) 

C 

NREC3-NEXREC 

INDEX ( 1 , DSNUM ) -NREC3 

WRITE ( 3 1 NREC3 )DSNAM, DATE. CONFIG, NMN, HAXNAA. NDV, NPAR, 

* ( ( PARMS ( I , J ) , J-l , 3 ) . 1-1 , NPAR ) . ( PARVAL (I ) . 1-1 . NPAR) 

C 

INDEX ( 2 , DSNUM ) -NREC3 

WRITE (3'NREC3) ( VAR(1 , J ) , J-l . 3 ) ,XLIM 
C 

READ( 1,1005) ( (A(M, I) , 1-1 , MAXNAA) ,M-1 , NMN) 

C 

DO 30 M-l.NMN 

30 WRITE ( 3 ' NREC3 ) ( A( M, I ), 1-1 , MAXNAA) 

C 

INDEX ( 3 , DSNUM ) -NREC3 

WRITE ( 3 ' NREC3 ) (VAR(2 , J ) , J*1 , 3 ) , APLIM 
C 

READ(1, 1005) ( (A(M, I) ,1-1, MAXNAA) , M-l.NMN) 

C 

DO 40 M-l , NMN 

40 WRITE ( 3 ' NREC 3 ) NAAX(M) , (A( M. I ). 1-1 , MAXNAA) 

C 

READ( 1,1005) ( DUM( I ) , 1=1 , NMN ) 

C 

I NDEX ( 4 , DSNUM ) -NREC3 

WRITE ( 3 ’ NREC 3 ) ( DUM ( I ) , I— 1 , NMN ) , ( ( RUNNO ( I , J ) , 1-1 . NMN ) , J-l , 2 ) 

C 

DO 60 K-l , NDV 
KK-K+4 

INDEX ( KK , DSNUM ) -NREC 3 

WRITE ( 3 ' NREC3 ) vDVAR(K, I ) , 1-1 , 3 ) , (DLIM(K, I ) . 1-1 , 2) 

C 

READ(1, 1005) ( (A(M, I) ,1*1, MAXNAA), M-l, NMN), (DUM(I).I-l, NMN) 

C 

DO 50 M-l.NMN 

50 WRITE(3'NREC3)(AfM, I). 1-1, MAXNAA) 

60 CONTINUE 

NEXREC-NREC3 
REWIND 1 

SUMMARY? 

65 WRITE(5, 5001 ) SETID 
READ (5,5002) ANSWER 
IF ( ANSWER. EQ. YES) CALL SUMMRY 

GET NAMES OF INDEPENDENT VARIABLES 

DO 68 1-1,2 

NREC3-INDEX (I +1 , DSNUM ) 

READ (3 ’NREC3) ( VAR (I , J ) , J-l , 3 ) 

GO CONTINUE 
C 

RETURN 
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70 CONTINUE 
FLAG-1 
REWIND 1 
RETURN 

FORMATS********************* 

1001 FORMAT (6X.3A2) 

1002 FORMAT(6X, 33A2, 12X, 416) 

1003 FORMAT (66A2) 

1004 FORMAT (2216) 

1005 FORMAT ( 8E 15.8) 

5001 FORMAT (/26H$DO YOU WANT A SUMMARY OF , 3A2,2H? ) 

5002 FORMAT (Al) 


END 


oono 


SUBROUTINE FINDER 


THIS SUBROUTINE GETS BLOCKS OF DATA FROM 
A DATABASE FOR USE BY SUBROUTINE ABC 

IMPLICIT INTEGER (A-Z) 

COMMON /IDENTS/ SETID( 3 ) , VAR(2 , 3 ) , IVID(3 ) , DVID(3 ) , ACPARM( 5 ) , 

* NAMES (10, 3 ) 

COMMON /COUNT/ MODE, NCURV, NSETS, NAPARM, NFEE, MODP, NPT, FLAG, 

* IVFLAG, ACFLAG, MACP 
COMMON /FILE3/ NREC3 , NEXREC 

COMMON /NDEX/ INDEX(14, 100) ,DSLIST(100, 3) ,DSNUM, SETCNT 
COMMON /BETA/ A( 10, 20 > , B( 10. 20 ) ,C( 10, 20 ) ,NBP( 10 ) . 

* SETNAM( 3 ) .CONFIG (24) , DVNAM( 3 ) , NMN, MAXNAA, NPARM, 

* NAMPAR (3,10), PVAL (10), FEEVAL, TMACH, TALPH 

C 

REAL A , B , C , PVAL , FEEVAL , TMACH , TALPH , ACPARM 
C 

DIMENSION FEE ( 3 ) , DATE ( 6 ) 

C 

DATA FEE /2HPH,2HI , 2H / 

C 

NREC 3=INDEX ( 1 , DSNUM) 

READ ( 3 ' NREC3 ) SETNAM, DATE, CONFIG, NMN, MAXNAA, NDV, NPARM. 

* ( (NAMPAR(I, J) ,1=1,3) ,J-1, NPARM) , ( PVAL(J ) ,J*1 , NPARM) 
C 

IF ( MODE . NE . 3 ) GO TO 200 
C 

M=0 

DO 150 N=l, NPARM 

IF ( NAMPAR ( 1 , N ) . NE . FEE ( 1 ) ) GO TO 50 
IF ( NAMPAR ( 2 , N ) . NF, . FEE ( 2 ) ) GO TO 50 
IF ( NAMPAR ( 3 , N ) . NE . FEE ( 3 ) ) GO TO 50 
C 

FEEVAL=PVAL(N) 

GO TO 150 
50 M*M+1 
C 

DO 100 1*1,3 

NAMPAR ( I , M ) =NAMPAR ( I , N ) 

100 CONTINUE 
C 

150 CONTINUE 
C 

N D ARM=M 

C 

200 CONTINUE 

NPEC3-INDEX ( 2 , DSNUM ) +1 
DO 250 1*1, NMN 

READ(3 ’NPEC3) (A( I , J ), J*1 , MAXNAA) 

250 CONTINUE 
C 

NREC3-INDEX ( 3 , DSNUM ) +1 
DO 300 1*1 , NMN 

READ (3 ’NREC3) NBP ( I ) , ( B( I , J ) , J*1 , NBP ( I ) ) 

300 CONTINUE 
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c 

DO 350 K-l.NDV 
NREC3-INDEX(K+4, DSNUM) 

READ ( 3 ' NREC3 ) DVNAM 

IF (DVNAH(1 ) .NE.DVID(1 ) ) GO TO 350 
IF (DVNAM( 2 ) .ME. DVID(2 ) ) GO TO 350 
IF (DVNAM( 3 ) . EQ. DVID( 3 ) ) GO TO 400 
350 CONTINUE 

NOT FOUND - SET FLAG AND RETURN 

FLAG-1 
RETURN 
C 

400 CONTINUE 

DO 450 I-l.NMN 

READ ( 3 ' NREC3 ) ( C( I , J ) , J-l , MAXNAA) 

450 CONTINUE 
C 

RETURN 

END 
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SUBROUTINE MYSTRO 

THIS SUBROUTINE PROMPTS THE USER FOR PLOT 
SPECIFICATIONS AND CONTROLS THE PRODUCTION OF 
THE CURRENT FRAME. 

IMPLICIT INTEGER (A-Z) 

COMMON /COUNT/ MODE , NCURV, NSETS, NAPARM, NFEE, MODP, NPT, FLAG, 

* I VFLAG , ACFLAG , MACP 

COMMON /IDENTS/ SETID( 3 ) , VAR (2 , 3 ) , IVID( 3 ) , DVID( 3 ) , ACPARM( 5 ) , 

* NAMES (10,3) 

COMMON /NDEX/ INDEX ( 14 , 100 ) , DSLIST ( 100, 3 ) , DSNUM , SETCNT 
COMMON / CURVS/ X ( 20, 5 ) , Y( 20, 5 ) , NP ( 5 ) , AVAL( 5 ) , CPAR 
COMMON /TEXT/ IDCURV (6,5), IDSETS (3,10), IDCFIG (24,5), 

* IDP\RM( 3 ,10,5) ,NPARMS(5 ) , PARVAL(10, 5 ) 

COMMON /BETA/ A( 10, 20 ) , B( 10, 20 ) , C( 10, 20 ) , NBP (10 ) , 

* SETNAM( 3 ) , CONFIG (24 ) , DVNAM( 3 ) , NMN, MAXNAA, NPARM, 

* NAMPAR(3 , 10 ) , PVAL(10 ) , FEEVAL, TMACH, TALPH 

DIMENSION FELIST(IO) 

REAL ACPARM, AVAL, CPAR, X, Y, A , B, C, PARVAL, PVAL, FEEVAL, TMACH, TALPH 

DATA YES / 1HY/ , 

* MOCK / 2HMA/ , 

* ALFA / 2HAL/ , 

* CRT / 2HCR/ , 

* FEE / 2HPH/ , 

* ZERO/O/ 

IF (IVID(1 ) .EQ.FEE) GO TO 300 

I.V. IS ALPHA, MACH NUMBER, OR CRT 
PROMPT FOR DATASET NAME 

50 NSETS*NSETS+1 
WRITE (5, 5001 ) 

READ (5,5002) SETID 

CALL CREATO 

IF ( FLAG. EQ. ZERO) GO TO 100 
FLAG=0 

WRITE (5,5003) SETID 
GO TO 50 

PROMPT FOR DEP. VAR. 

100 WRITE(5, 5005) 

READ (5,5002) DVID 


CALL FINDER 

IF ( FLAG. EQ. ZERO) GO TO 250 
FLAG-0 

WRITE (5,5018) DVID, SETID 
GO TO 100 
250 CONTINUE 
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IF ( IVID( 1 ) . EQ.VAR(1 , 1 ) ) GO TO 260 
IF ( IVID(1 ) . EQ. VAR(2 , 1 ) ) GO TO 270 

I.V. NOT FOUND - RETURN 

FLAG-1 

RETURN 

260 I VF LAG-1 
ACFLAG-2 
MODE-1 
GO TO 280 

270 I VF LAG-2 
ACFLAG-1 
MODE-2 

280 CONTINUE 

WRITE ( 5 , 5006 ) (VAR (ACFLAG, I ) , 1-1 , 3 ) 

READ(5 , 5011 ) NAPARM 

WRITE ( 5 , 501 6 ) NAPARM , ( VAR ( ACFLAG , I ) , 1-1 , 3 ) 

READ ( 5 , 5008 ) ( ACPARM ( I ) , 1-1 , NAPARM ) 

PROCESS DATA ACCORDING TO MODE 

CALL ABC 

ANYMORE ON THIS FRAME? 

WRITE (5, 5009) 

READ (5, 5004) ANSWER 

IF ( ANSWER. NE. YES) GO TO 290 

SAME DATASET? 

WRITE (5, 5015) 

READ ( 5 , 5004 ) ANSWER 
IF (ANSWER. EQ. YES) GO TO 100 
GO TO 50 

290 RETURN 

PHI IS I.V. - PROMPT FOR NUMBER OF DATASETS 

300 MODE-3 

WRITE(5, 5010) 

READ (5,5011) NFEE 

READ DATASET NAMES 

WRITE (5, 5012) NFEE 
DO 350 N»1 , NFEE 

350 READ (5,5002) ( NAMES (N, I ) , 1-1 , 3) 
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DO 380 N-l , NFEE 
DO 360 1-1,3 

SETID ( I )-NAMES (N , I ) 

360 CONTINUE 
C 

CALL CREATO 
FELIST (N ) -DSNUM 
380 CONTINUE 
C 

C READ DEP. VAR. NAME 

C 

WRITE (5, 5005) 

READ (5,5002) DVID 
C 

C READ INVARIANT PARAMETER NAME AND VALUE 

C 

400 MODP-O 

WRITE (5, 5013) ( (VAR(I , J ) , J-l , 3 ) , 1-1 , 2 ) 

READ (5,5002) CONPAR 
IF ( CONPAR . EQ . VAR (1,1)) MODP-1 
IF (CONPAR. EQ.VAR(2,1)) MODP-2 
C 

C READ ACTIVE PARAMETRIC VALUES 

C 

450 WRITE (5, 5014) 

IF (MODP.EQ.l) GO TO 500 
READ (5,5008) TALPH 
CPAR-TALPH 

WRITE(5, 5006) ( VAR( 1 , J ) , J-l , 3 ) 

READ (5,5011) NAPARM 

WRITE (5,5016) NAPARM , ( VAR( 1 , J ) , J-l , 3 ) 

GO TO 550 

500 READ (5,5008) TMACH 
CPAR-TMACH 

WRITE(5, 5006 ) (VAR(2, J) , J-l , 3) 

READ (5,5011) NAPARM 

WRITE (5,5016) NAPARM, (VAR (2, J) , J-l , 3 ) 

550 READ (5,5008) ( ACPARM( I ), 1-1 , NAPARM) 

C 

DO 650 N-l, NFEE 
NSETS-NSETS+1 
C 

DSNUM-FELI ST ( N ) 

C 

CALL FINDER 
C 

NPT-N 
CALL ABC 
650 CONTINUE 
C 

• NCURV-NAPARM 
CALL SORTOR 
RETURN 
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FORMATS 


5001 FORMAT (/28H$ENTER DATASET IDENTIFIER - ) 

5002 FORMAT (3A2) 

5003 FORMAT ( /IX, 29HDATAMANIA CAN'T FIND DATASET , 3A2/) 

5004 FORMAT (Al) 

5005 FORMAT (/33H$ENTER DEPENDENT VARIABLE NAME - ) 

5006 FORMAT (/27H$ENTER NUMBER OF VALUES OF , 3A2.17H TO BE PLOTTED - ) 

5008 FORMAT (5F10.0) 

5009 FORMAT ( / 3 OH $ ANY MORE DATA ON THIS FRAME? ) 

5010 FORMAT (/44H$ENTER NUMBER OF PHI VALUES TO BE PLOTTED - ) 

5011 FORMAT (II) 

5012 FORMAT (/ IX, 6HENTER ,I1,14H DATASET NAMES/ 

* IX, 32HEACH FOLLOWED BY CARRIAGE RETURN) 

5013 FORMAT (/38H$ENTER PARAMETER TO BE HELD CONSTANT (,3A2,4H OR , 

* 3A2 , 4H ) - ) 

5014 FORMAT (/37H$ENTER VALUE OF CONSTANT PARAMETER - ) 

5015 FORMAT (/19H$USE SAME DATASET? ) 

5016 FORMAT (/6H5ENTER, 12, 11H VALUES OF ,3A2, 

* 26H (SEPARATE WITH COMMAS) - ) 

5018 FORMAT ( /IX, 3A2 , 28H IS NOT A DEPENDENT VARIABLE 

* /IX , 11HIN DATASET ,3A2,13H. TRY AGAIN./) 


END 
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SUBROUTINE NSTRUC 

THIS SUBROUTINE WRITES INSTRUCTIONS FOR THE USE OF PROGRAM 
DATAMANIA 

WRITE (5, 5001) 

WRITE (5, 5002) 

WRITE (5, 5003) 

5001 FORMAT ( ///IX , 48H DATAMANIA IS DESIGNED TO READ DATAMAN TAPES 
/1X.55HAND PRODUCE PLOTS OF COEFFICIENT DATA AS A FUNCTION OF 
/IX, 52HEITHER ANGLE OF ATTACK (ALPHA), MACH NUMBER, OR ROLL 
/IX, 55HANGLE (PHI). IN ORDER TO DO THIS, YOU NEED TO KNOW THE 
/IX, 54HNAMES OF THE DATASETS WHICH CONTAIN THE INFORMATION TO 
/IX, 11HBE PLOTTED.) 

5002 FORMAT ( /IX, 51H AT THE BEGINNING OF EACH FRAME, DATAMANIA WILL 

* /IX, 55HPROMPT FOR THE INDEPENDENT VARIABLE NAME. YOUR RESPONSE 

* /IX, 49HOF ALPHA, MACH, OR PHI WILL DETERMINE THE MODE OF 

* /IX, 5 5H PROGRAM OPERATION. FOLLOWING THIS YOU WILL BE PROMPTED 

* /IX, 45HFOR INFORMATION NEEDED TO COMPLETE THE FRAME.) 

5003 FORMAT ( /IX, 49H AS EACH FRAME IS COMPLETED YOU WILL BE ASKED 

* /IX, 52HFIRST IF YOU WANT A TABULATION OF THE CURRENT FRAME, 

* /IX, 50HAND THEN, IF YOU WANT TO CREATE A PLOT FILE OF THE 

* /IX, 6HFRAME. / ///) 

RETURN 

END 


SUBROUTINE PLOTOR 

THIS SUBROUTINE CREATES AN SGP PLOT FILE 
FROM THE DATA PRODUCED BY SUBROUTINE ABC 

IMPLICIT INTEGER (A-Z) 

C 

COMMON / CURVS/ X ( 20, 5 ) , Y( 20 , 5 ) , NP (5 ) , AVAL( 5 ) , CPAR 
COMMON /TEXT/ IDCURV(6 , 5 ) , IDSETS ( 3 , 10 ) , IDCFIG (24, 5 ) , 

* IDPARM ( 3,10,5), NPARMS ( 5 ) , PARVAL( 10, 5 ) 

COMMON /COUNT/ MODE, NCURV, NSETS, NAPARM, NFEE, MODP, NPT, FLAG, 

* IVFLAG, ACFLAG.MACP 

COMMON /IDLNTS/ SETID ( 3 ) , VAL( 2 , 3 ) , IVID( 3 ) , DVID( 3 ) , ACPARM( 5 ) , 

* NAMES (10, 3) 

COMMON /BSCALE/ XMAX, XMIN, YMAX, YMIN 

C 

DIMENSION CNTRL(2, 5) ,XLABEL(4) ,YLABEL(4) ,SYMB(5) 

DIMENSION NAME (3 ) , SHAP1 ( 12 ) , SHAP2 (12 ) , PARM( 3 ) , CURV1 (3 ) , CURV2 ( 3 ) 
DIMENSION ZLABEL(4 ) , PLABEL(4 ) 

DIMENSION WLABEL ( 4 ) , SHAPE (24) 

C 

REAL X.Y.PARVAL, AVAL, CPAR, XMAX, XMIN, YMAX, YMIN 
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c 


c 

c 

c 


c 


c 

c 

c 


c 

c 

c 


c 

c 


c 


c 

c 

c 

c 


DATA CNTRL /350, 0, 180, 340, 90, 260, 24, 213 , 10, 170/ 
DATA NDP , I PUT , JPNT , STAR /20,2,1,2H* / 

DATA YLABEL /2H *,2H , 2H , 2H / 

DATA PLABEL /2H , 2HPH.2HI , 2H* / 

DATA SYMB /2H0 , 2H1 , 2H2 , 2H3 , 2H4 / 

DATA YES /1HY/ 

SET SCALING OPTIONS 

IOPT-2 

WRITE (5, 5001) 

READ (5,5002) ANSWER 

IF ( ANSWER. EQ. YES) GO TO 9* 

IOPT-1 

WRITE (5, 5003) 

READ (5, 5002) ANSWER 
IF (ANSWER. EQ. YES) Go TO 99 

READ USER DEFINED SCALING VALUES 

WRITE (5, 5004) 

READ( 5 , 5005 ) XMIN.XMAX 
WRITE (5, 5006) 

READ (5, 5005) YMIN.YMAX 

SET LABELS 

99 CONTINUE 

IF ( MODE . GT . 2 ) GO TO 4 
DO 1 1-1,3 

XLABEL( I )-VAR( IVFLAG, I ) 

Z LABEL ( I ) -VAR ( AC FLAG , I ) 

1 CONTINUE 

XLABEL ( 4 ) -STAR 
ZLABEL(4)-STAR 
GO TO 7 


4 CONTINUE 

DO 5 1-1,4 

XLABEL ( I ) -PLABEL ( I ) 

5 CONTINUE 

DO 6 1-1 , NCURV 

NP( I )-NFEE 

6 CONTINUE 

7 CONTINUE 

CALL PLOTI (X , Y, IOPT, NP, NDP, NCURV, IPNT, JPNT, XLABEL, YLABEL) 

IF ( MODE . EQ . 3 ) GO TO 90 

XINC-CNTRL( 2 , NSETS ) 

XLOC -CNTRL ( 1 , NSETS ) 
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c 

c 


c 


c 


c 


c 


c 


c 


DO 50 N-l , NSETS 
DO 10 1-1,3 

10 NAME(I)-IDSETS(I,N) 

DO 20 1-1,12 

SHAP1 ( I ) -I DCF 10 ( I ,N) 

20 SHAP2(I )-IDCFIG(I+12,N) 

CALL TXT(XLOC, 1100, NAME, 6,0) 

CALL TXT(XLOC, 1085, SHAP1, 24,0) 

CALL TXT(XLOC, 1070, SHAP2, 24,0) 

YLOC-1055 

NUMLOC-XLOC+50 

DO 40 M— 1 , NPARMS ( N ) 

DO 30 1-1,3 

30 PARM ( I ) - I PPARM ( I , M , N ) 

CALL TXT (XLOC, YLOC, PARM, 6,0) 

CALL NUMB ( NUMLOC , YLOC , PARVAL (M,N),6H(F6.2),0) 

YLOC-YLOC-1 5 
40 CONTINUE 

XLOC-XLOC+XINC 


C 

50 CONTINUE 
C 

MAX-0 

DO 60 I-l.NCURV 

MAX— MAXO ( MAX , NPARMS ( I ) ) 

60 CONTINUE 
C 

XLOC-200 

YLOC-1 040-MAX* 15 
C 

CALL TXT(XLOC, YLOC, 28HSYMBOL DV DATASET , 28 , 0 ) 

CALL TXT (300, YLOC, ZLABEL, 6,0) 

YLOC -YLOC -8 

CALL TXT(XLOC, YLOC, 28H — ,28,0) 

C 

DO 80 N-l , NCURV 
C 

DO 70 1-1,3 

CURV1 ( I )-IDCURV ( I , N ) 

CURV2 ( I )-IDCURV( 1+3 , N ) 

70 CONTINUE 
C 

YLOC-YLOC-1 5 

CALL TXT (215, YLOC , SYMB ( N ) ,2,0) 

CALL TXT(265, YLOC, CURV1, 6,0) 

CALL NUMB ( 300, YLOC, AVAL (N ) , 6H ( F6 . 2 ) , 0 ) 

CALL TXT (360, YLOC, CURV2, 6,0) 

80 CONTINUE 
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GO TO 450 
C 

C IND. VAR. IS PHI 

C 

90 CONTINUE 
C 

c 

DO 92 1-1,3 

ZLABELU )-VAR(MACP, I ) 

WLABEL ( I ) -VAR ( MOD? , I ) 

92 CONTINUE 
C 
C 
C 

DO 100 1-1,3 

CURV1 (I )— IDCURV(I, 1 ) 

CURV2 ( I )«IDCURV( 1+3 , 1 ) 

100 CONTINUE 
C 

CALL TXT (350, 1100, CURV1, 6,0) 

CALL TXT (400, 1100, 5H VS. ,5,0) 

CALL TXT (440, 1100, PLABEL,6,0) 

C 

DO 150 1-1,24 

SHAPE ( I )— IDCFIG (1,1) 

150 CONTINUE 
C 

CALL TXT (200, 1070, If .CONFIGURATION - ,16,0) 
CALL TXT (320,1070.' IAPE. 48,0) 

C 

CALL TXT (125, 1030, 6HSYMBOL, 6, 0) 

CALL TXT (175, 1030, ZLABEL,6,0) 

CALL TXT (300, 1030, 17HPARAMETRIC VALUES, 17,0) 
CALL TXT (600, 1030, 8HDATASETS, 8,0) 

C 

YLOC-1010 

C 

DO 200 N— 1 , NCURV 

CALL TXT (140, YLOC, SYMB(N) ,2,0) 

CALL NUMB (180, YLOC , AVAL (N),6H(F5.2),0) 
YLOC-YLOC-1 5 
200 CONTINUE 
C 

YLOC-1010 

C 

DO 300 N— 1 , NPARMS ( 1 ) 

DO 250 1-1,3 

PARM( I )— IDPARM( I , N, 1 ) 

250 CONTINUE 
C 

CALL TXT (310,YLOC,PARM,6,0) 

CALL NUMB ( 3 6 0 , YLOC , PARVAL ( N , 1 ) ,6H(F6.2) ,0) 
YLOC-YLOC-1 5 
300 CONTINUE 
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c 

YLOC-1010 

C 

DO 400 N-l , NSETS 
C 

DO 350 1-1,3 

NAME (I )— IDSETS(I ,N) 

350 CONTINUE 
C 

CALL TXT (610, YLOC, NAME, 6.0) 

YLOC-YLOC-15 
400 CONTINUE 
C 

CALL TXT ( 400 , 800 , WLABEL, 8,0) 

CALL NUMB(445,800,CPAR,6H(F5.2) .0) 
C 

450 CALL ENDFRM 
C 

RETURN 


FORMATS********************* 

5001 FORMAT ■( /22H$USE DEFAULT SCALING? ) 

5002 FORMAT (Al) 

5003 FORMAT ( / 30H$USE SAME SCALE AS LAST PLOT? ) 

5004 FORMAT (/ 36H ENTER MINIMUM AND MAXIMUM VALUES OF 

* /45H$INDEPENDENT VARIABLE (SEPARATED BY COMMA) - ) 

5005 FORMAT (2F1 0.0) 

5006 FORMAT (/36H ENTER MINIMUM AND MAXIMUM VALUES OF 

* /43H$DEPENDENT VARIABLE (SEPARATED BY COMMA) - ) 

C 

£******************************************** 

C 

END 
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SUBROUTINE SORTOR 

C 

COMMON /CUR VS/ X(20. 5) ,Y(20, 5) ,NP(5 ) , AVAL(5) ,CPAR 
COMMON /COUNT/ MODE, NCURV, NSETS, NAPARM, NFEE, MODP.NPT, SAMSET 
C 

NPM-NFEE-1 

C 

DO 300 I-l.NPM 
JJ-I+1 

DO 200 J«J J , NFEE 

IF (X(I, 1) .LT.X(J,1) ) GO TO 200 

TEMP*X(J,1) 

X(J.1)=X(I,1) 

X ( 1 , 1 )=TEMP 

DO 100 K=1 , NCURV 

TEMP*Y(J,K) 

Y(J,K)=Y(I,K) 

Y(I,K)-TEMP 
100 CONTINUE 

200 CONTINUE 
300 CONTINUE 
C 

DO 400 J*2 , NCURV 

DO 400 I®1 , NFEE 

X(I,J)-X(1,1) 

400 CONTINUE 
C 

RETURN 

END 


SUBROUTINE SUMMRY 

THIS ROUTINE PRODUCES INFORMATION ABOUT THE DATASET 
CURRENTLY BEING SEARCHED 

IMPLICIT INTEGER (A-Z) 

C 

COMMON /FILE3/ NREC3 , NEXREC 

COMMON /NDEX/ INDEX (14, 100) ,DSLIST(100, 3) , DSNUM , SETCNT 
C 

DIMENSION NAME ( 3 ) .DATE (6) .BODY (24) , PARMS( 10, 3 ) , DV(3 ) 

C 

REAL PARVAL (10), MOX ( 20 ) , MIN , MAX 
C 

NREC3 r j uDEX ( 1 , DSNUM ) 

READ (3 'NREC3 ) NAME, DATE, BODY, NMN, MAXNAA, NDV.NPAR, 

* ( (PARMS( I, J ) , J*1 , 3 ) , 1=1 ,NPAR) , ( PARVAL ( I ) , 1*1 , NPAR) 

WRITE (5, 5001) NAME, DATE, BODY 
C 

DO 10 1*1 , NPAR 

10 WRITE(5, 5002 ) ( PARMS( I , J ) , J*1 , 3 ) , PARVAL ( I ) 
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c 

WRITE (5, 5003) 

C 

NREC3«INDEX(2, DSNUM)+1 
C 

DO 30 I-l.NMN 

READ (3 ’NREC3) (MOX(L) , L»1 , MAXNAA ) 

MIN“MOX( 1 ) 

MAX-MOX(l) 

C 

DO 20 J-2 , MAXNAA 

MIN“AMINl (MIN, MOX( J ) ) 

MAX=AMAX1 (MAX. MOX( J ) ) 

20 CONTINUE 
C 

WRITE (5, 5004) I, Mil, MAX 
30 CONTINUE 
C 

WRITE (5, 5005) 

C 

DO 40 I»1,NDV 
NREC3*INDEX(I+4, DSNUM) 

READ ( 3 'NREC3 ) DV 
WRITE (5, 5006) DV 
40 CONTINUE 
C 

WRITE (5, 5007) 

C 

RETURN 

C 

£************************pO£MArpg*************** ********* 

C 

5001 FORMAT (//13X.11HSUMMARY OF , 3A2/1X, 45 ( 1H*) /IX, 1H* , 

* 32X, 8HCREATED , 6A2/1X, 1H* , 5X, 16HCONFIGURATION - , 

* 24A2/1X, 1H*/1X, 1H*,14X, lOHPARAMETERS/lX, 1H*. 

* 13X, I2(1H-) ) 

5002 FORMAT (IX, 1H*, 10X, 3A2, 3H = ,F6.2) 

5003 F0RMAT(1X, 1H*/1X, 1H*, 13X, 12HMACH NUMBERS/1X, 1H* , 13X, 

* 12(lH-)/lX, 1H*,6X, 5H INDEX, 10X, 5H RANGE ) 

5004 FORMAT ( IX, 1H*,6X, 12, 6X, 5HFROM ,F6.3,4H TO ,F6.3) 

5005 FORMAT (IX, 1H*/1X, 1H*,10X, 19HDEPENDENT VARIABLES/lX, 1H*,9X, 21(1H-) ) 

5006 F0RMAT(1X, 1H*, 17X, 3A2) 

5007 FORMAT (IX, 1H*/1X, 1H*/1X, 45 ( 1H* ) ) 

******************************************************** 


END 
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ooooo no noon 


SUBROUTINE TABULA 
C 

IMPLICIT INTEGER (A-Z) 

C 

COMMON / CURVS/ X (20, 5 ) ,Y (20. 5 ) ,NP( 5 ) , AVAL(5 ) ,CPAR 

COMMON /TEXT/ IDCURV(6, 5) , IDSETS(3, 10) , IDCFIG('4, 5) . 

* IDPARM( 3.10,5), NPARMS ( 5 ) , PARVAL( 10. 5 ) 

COMMON /COUNT/ MODE , NCURV , NSETS , NAPARM , NFEE , MODP , NPT , FLAG, 

* I VFLAG ACFLAG # MACP 

COMMON /IDENTS/ SETID(3) ■ VAR(2, 3) , IVID(3) ,DVID(3 ) , ACPARM(5) , 

* NAMES (10, 3) 

C 

REAL X, Y , AVAL, CPAR , PARVAL 
C 

DIMENSION FEE ( 3 ) 

C 

DATA FEE /2H P,2HHI,2H / 

C 

WRITE (5, 5001) 

C 

IF (MODE.EQ.3) GO TO 450 


DO 300 N-l, NSETS 

WRITE ( 5 , 5003 ) ( IDSETS ( I . N ) , 1-1 , 3 ) , ( IDCFIG ( I , N ) , 1-1 , 24 ) 

WRITE (5, 5004) 

DO 250 M-l, NPARMS (N) 

WRITE ( 5 , 5005 ) ( 1DPARMS ( I , M , N ) , 1-1 , 3 ) , PARVAL (M,N ) 

250 CONTINUE 

300 CONTINUE 

WRITE (5, 5008) 

DO 400 N«l, NCURV 

WRITE ( 5 , 5006 ) ( VAR ( I VFLAG , J ) , J-l , 3 ) , ( I DCURV ( I , N ) , 1-1 , 3 ) , 

* (VAR(ACFLAG, J) ,J-1,3) ,AVAL(N), (IDCURV(I.N) ,1-4,6) 

DO 350 M»1,NP(N) 

WRITE(5, 5007 ) X(M,N) ,Y(M,N) 

350 CONTINUE 

400 CONTINUE 

WRITE (5, 5008) 

RETURN 


450 CONTINUE 

WRITE (5, 5011) 

DO 500 N-l, NSETS 
WRITE ( 5 , 5002 ) ( IDSETS (I ,N) , 1-1 , 3 ' 
500 CONTINUE 
C 

WRITE (5, 5012) (IDCFIG (1,1) ,1-1,24) 
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noon 


c 


MACP-2 

IF (M0DP.EQ.2) MACP-1 


WRITE (5, 5004) 

DO 750 M-l.NPARMS(l) 

WRITE ( 5 , 5005 ) ( I DP ARMS ( I . M, 1 ) , 1-1 . 3 ) , PARVAL ( M. I ) 

750 CONTINUE 
C 

WRITE (5, 5008) 

WRITE (5,5009) (VAR ( MODP , J ) , J-l , 3 ) , CPAR 
DO 850 N=1 , NCURV 

WRITE ( 5 , 5010 ) FEE, ( IDCURV (I,N) , 1*1, 3) , ( VAR (MACP , J ) , J«1 , 3 ) , AVAL (N ) 
DO 800 M*1 , NFEE 

WRITE(5, 5007) X(M, N) , Y(M, N) 

800 CONTINUE 
850 CONTINUE 
C 

WRITE (5, 5008) 

C 

RETURN 

C 

C* *****************pQjyiaij>g*** ************** ****** 

C 

5001 FORMAT ( / / / /IX, 16 ( 1H$ ) , 28H TABULATION OF CURRENT FRAME ) 

5002 F0RMAT(1X, 1H$, 16X, 3A2 ) 

5003 F0RMAT(/1X, 1H$/1X, 1H$/1X, 1H$, 4X, 10HDATASET - ,3A2/1X,1H$, 

* 4X, 16HCONFIGURATION - ,24A2) 

5004 FORMAT (IX, 1H$, 17X, 17HPARAMETRIC VALUES/ 

* IX, 1H$ , 17X, 17H ) 

5005 FORMAT (IX, 1H$, 17X, 3A2,4X,F6,2) 

5006 F0RMAT(/1X, 1H$/1X, lH$/lX, 1H$,6X, 3A2,8X, 3A2, 2X.4HFOR ,3A2, 

* 1H=,F6.2,6H FROM ,3A2) 

5007 FORMAT (IX, 1H$ , 2X, F10. 4, 3X.F10.4) 

5008 FORMAT (IX, 1H$/1X, 43 (1H$ ) ) 

5009 FORMAT ( /20X, 3A2 , 1H*, F5. 2/ ) 

5010 F0RMAT(1X, 1H$/1X, 1H$/1X, 1H$, 6X, 3A2, 8X, 3A2, 2X, 4HFOR ,3A2, 

* 1H*,F6.2) 

5011 FORMAT (IX, 1H$/1X, 1H$, 15H USING DATASETS ) 

5012 FORMAT (2H $/19H $ CONFIGURATION - , 24A2/2H $) 

C 

Q + ********************************************* 

C 

END 
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TABLE 1. MISSILE TESTS INCLUDED IN DATA BASE 


TEST NUMBER 

CALSPAN T17-093 
CALS PAN T17-123 
ARC 66-036 
AEDC V41A-C1A 
AEDC SF172/TF360 
AEDC TF-416 
ARC T183-11 


TITLE REFERENCE 

Canard Control 1 

Plume Effects 2 

Canard Control 3 and 8 

Canard Control 4 

Plume Effects 5 

Plume Effects 6 

High-Angle-of-Attack 7 
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TABLE 2. FORMAT AND NOMENCLATURE FOR DAT/ ’'IAN TAPES 


(a) Format Statements 


The DATAMAN tapes are written using these FORTRAN 
statements. 


WRITE (KTAPE, 100) FDSET , IDENT, DATE (1) , DATE (2) , (CONFIG(I), 

I = 1,8), ITYPE , MODE, NMN, MAXNAA, NDV, 
NPAR, NDM 


100 FORMAT (12A6, 716) 

WRITE (KTAPE, 101) (IDPVAR(I) , 1=1,2), (NDPVAR(I) , 1 = 1,2), 

(NAMPAR(I) , 1 = 1, NPAR), (FILEREF (I) , 

1 = 1,3), (REF (1,2) , 1 = 1,7), (NAMDM(I), 

I = 1,20) 

101 FORMAT (22A6) 

WRITE (KTAPE, 102) ( (IRUNNO (I, J) , I = 1,NMN) , J = 1,2), 

(NAAX(M) , M = 1,NMN) 

102 FORMAT (2216) 

WRITE (KTAPE, 103) ANGMIN, ANGMAX, (XMLIM(I), I = 1,2), 

(APLIM(I), 1 = 1,2), ( (DLIMIT(I, J) , 1 = 1, 
NDV), J » 1,2), (PARVAL(I), 1 = 1, NPAR), 

( REF (1,1) , I = 1,7), ( (XMACH(M,I) , 1=1, 
MAXNAA), M = 1 , NMN ) , ( (ALPHA (I ,M) , 1=1, 
MAXNAA), M = 1, NMN), (RNL(M) , M = 1,NMN) , 
( (DM (M, I ) , M = 1, NMN), I = 1,20) 

DO 10 K = 1, NDV 

10 WRITE (KTAPE, 103) ((D(I,M,K), I = 1, MAXNAA), M = 1, NMN), 

(GD (M,K) , M = 1, NMN) 


103 FORMAT (8E15. 8) 


TABLE 2. Continued 


(b) Description of DATAMAN Variable Names 


Variable 

Type 

Dimension 

Description 

FDSET 

A 

none 

Contains the characters *FDSET 

I DENT 

A 

none 

Six character label used to identify 
the dataset 

DATE 

A 

2 

Dataset creation date 

CONFIG 

A 

8 

Configuration of model tested 

ITYPE 

I 

none 

Data type (1-raw data) and 
2 - nominalized data) . Has value 
1 for all nine tapes. The dependent 
variable data are actually presented 
in coefficient form. 

MODE 

I 

none 

Data mode (1- pitch angle varies) 
and 2 - yaw angle varies) . Has value 
1 for all nine tapes. 

NMN 

I 

none 

Number of Mach numbers tested (£20) 

MAXNAA 

I 

none 

Maximum number of alphas per Mach 
number (£50) 

NDV 

I 

none 

Number of dependent variables (£.10) 

NPAR 

I 

none 

Number of parameters (£10) . The 
parametric values are constant for 
each dataset. 

NDM 

I 

none 

Number of univariate dependent 
variables (£20) . Has value 0 for 
all nine tapes. 

I DP VAR 

A 

2 

Names of independent variables. 
Contains the characters MACH and 
ALPHA for all nine tapes. 

NDPVAR 

A 

10 

Names of dependent variables 

NAMPAR 

A 

10 

Names of parameters 
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TABLE 2. Concluded 


Variable Type Dimension Description 


FILREF 

A 

3 

Document references. Blank for all 
nine tapes. 

REF 

F,A 

(7,2) 

Reference data values and units. 

NAMDM 

A 

20 

Names of univariate dependent 
variables. Blank for all nine tapes. 

IRUNNO 

I 

(2,10) 

Run and rerun numbers associated with 
each Mach number. 

NAAX 

I 

20 

Number of alphas associated with each 
Mach number. 

ANGMIN 

F 

none 

Lower limit used in gradient compu- 
tation. Set to zero on all nine 
tapes. 

ANGMAX 

F 

none 

Upper limit used in gradient compu- 
tation. Set to zero on all nine 
tapes. 

XMLIM 

F 

2 

Lower and upper limits for Mach 
number. 

APLIM 

F 

2 

Lower and upper limits for alpha. 

DLIMIT 

F 

(10,2) 

Lower and upper limits for dependent 
variables. 

PARVAL 

F 

10 

Parametric values. 

XMACH 

F 

(20,50) 

Mach numbers. 

ALPHA 

F 

(50,20) 

Angles of attack. 

RNL 

F 

20 

Reynolds numbers. 

DM 

F 

(20,20) 

Univariate dependent variable data. 
Set to zero on all nine tapes. 

D 

F 

(20,50, 

10) 

Dependent variable data. 

GD 

F 

(20,10) 

Gradient dependent variable data. 
Set to zero on all nine tapes. 
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Figure X. * Structure of DATAMANIA random access file 



012440 01 OCT 74 COOT lUOUTIlO 4 (M1C1T1) Mil M-IU J 7 t t 

IET4 II II 12 14 II -1 12-4 NK MI-T 

-Lit ZliSi -I.H Lit -I.H >.W >48.11 


| Segment 1 




1.844 

JLilIL. 


1 .502 

jl m 


ttttti 

1.104 

JjSSI 


1. 140 
1.448 


T7355 

1.448 


TTS44 

1.448 


77154 

1.448 


^ Segment 2 


mTh* -2.174 12.824 1 

7 - 2 . 1 J 2 i.W — r.i*4 J.itS 4.111 ♦.4il tf.W 

7 -J.I74 -1.174 4.414 1.444 4.084 ».I13 12.210 


2.44 2.44 | 71 78 4 4 


} 


S«gm«nt 3 
8«9Mnt 4 


C* -0.4M 2.01 • I 

-0.444 -4.044 4.1 Ji 57731 0.171 77427 2.011 

-4.17? 4.447 4.114 4.403 0.122 1.114 1.471 


IS 057 1.445 I 

4.414 4.412 Till 4.4I> 0.4» 5.4» 6.445 

4.148 4.144 4.334 4.141 0.274 4.114 0.112 


£8 U3Z I 

4.144 -4.244 -0.112 -4.474 -1.114 -1.742 -2.124 

4.287 4.444 -4.013 -4.133 -4.427 -4.748 -1.401 

CT : OJI 4.447 | 

"-r.iSI -5.5JJ -4. 475 -5.5ii -4.4m -5.55J 057 

-4.024 -0.427 -4.423 -0.014 -0.411 -0.412 -0.018 


Tr5 hT553 0751 

4.044 57575 0.044 0.488 57537 5.51 t -5.554 

0.437 0.044 0.040 0.037 0.031 4.023 0.027 


Segment 5 


£18 -0.018 0.000 I 

•4.408 -0.044 -0.008 -0.007 -O.OtO -0.012 

0.000 -0.004 -0.003 -0.404 -0.407 -0.407 


•0.018 

-0.001 


C»M -0.004 0.004 I 

4.444 57553 0.004 4.001 57555 -5.55i - 4.004 

4.403 0.004 0.004 4.402 0.400 0.403 -4.003 


CTI -0.014 0,010 1 

-0.010 -0.407 -4.003-0.004 -0.403 4.404 -4.002 

-0.401 0.002 0.042 4.002 0.001 0.010 4.004 


cic r o7o73 5.124 i 

-0.074 -0.032 -0.017 0.004 57537 57517 57733' 

-0.047 -0.017 -0.004 0.012 0.434 0.042 0.003 


Figure 2. - Example of DATAMANIA random acceaa dataeet. 
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(g) Army F2 tail fins. 


Figure 3 Continued. 










f icient 


















OfiL T17-123 (BFP BOD*. 
FIN IN FORWARD P06IT:CTJ 
BETA e.00 

PHI 0.00 

MACH 1.25 


RXE013 

AEDl TF3B0 BODY FIN. BF 
1 

BETA 0.00 

PHI 0.00 

FINPOS 3.00 

MACH 1.25 


*8077 

on. T17-123 (t'D BOD', 
FIN IN FORWAF: P06IT IQf. 

BETA 0.00 

PHI 0,00 

MAO! 1.25 


S^BQl Qv cgL. Q9I9SET 

0 CLM 0.00 RX9038 

1 CLM e.01 *F013 

2 CLM 6.00 RXE013 

3 CLM £.00 RX9037 







Normal - Force Coefficient 








0L9PRN CflNPRD CCWiROL, 
B1N1C2T2 

BETA 0.00 

phioc e .00 
m\X9L 0.00 
OCPQS 3.00 
TflLPOS 1 00 
DOC1 0.60 
DOC2 -3.00 

0003 e .00 

DOOM -3 .00 


qj^000 

CfL$°fN CfiN^RO OOT3X, 
BIN1C2T2 
BETP 0.00 
PHIOC 0.00 
PHITPL 0.00 
0CE3S 3.00 

hlpos i.0e 
0001 0.00 
DCNP2 0.00 
0003 0.00 

OOtM 0.00 


wee: 

QPL9KW CPNP : 3 CQNTF'X, 

BETP 0 01 
PHICNj 0.01 
PHlTf_ 0.® 

0CPC3 3.00 
TflLPOS 1.00 

ooc: e.0r 

DCNK 6.00 
DCND3 0.0C 
DOOM 6.0C 


RUEB94 

CQLSPOP CflNPRP COK~^, 
BlNICJTc 
BE TO 0.00 

PHIOC 0.00 
PHITJV 0.09 
CNOPOS 3.00 
TPLPCS 1.00 
OOC1 0.00 
OOC2 lb. 00 
DOC3 0.00 
DOOM .S .00 


Dataset 

AUC 086 
AUE 088 
AUE 091 
AUE 094 


Data 


♦ 

A 


Missile 2 

o 

□ 

o 

A 



Angle of Attack, Degrees 
(a) ^ - 0; inline tail fins. 

Figure 8.- Comparison of pitching-moment coefficient data 
and predictions from Program MISSILE 2 for Burt’s 
Canard-Conf iguration tasts; M w ■ 0.9. 



ORIGINAL PAGE ! r 
r ^ry)P QUALlT 


flUEiei 

AUE102 


HUE 105 


AUE 107 


CPLP°flN C t SW CC TROL. 

ccLf^ aw control. 

CALSPPN VhPPQ C0N7F3L, 

CALSMN CfHfD CCK^OL, 

B :C2T2 

Bit :C2T2 

B1N1C2T2 

B1N1C2T2 

BE- 0.W 

BETA 

0.00 

BETA 

0.00 

BETA 

0.00 

PH: _ ‘0 45. 20 

PHICND 

45. Cd 

PHICM) 

M5.00 

PHIOO 

45.00 

ph:t3_ 020 

°HITQL 

0.00 

PHITOL 

0.00 

PHITRL 

0.00 

Of'OS 3 20 

OOPOS 

3.00 

CNDPOS 

3.00 

OCPOS 

3.00 

T ft -OS 1 20 

Tft-POS 

1.00 

TflLPQS 

1.00 

TALPOS 

1 .00 

Dc:i -3.00 

DODi 

0.00 

0CND1 

6.00 

DCND1 

1500 

DC 02 -3.00 

DC?v03 

0.00 

DO02 

6.00 

DOD2 

15.00 

DO .23 -3.00 

0003 

0.00 

0003 

6.00 

0003 

15.00 

DC 2)4 -3.00 

DCND4 

0.00 

DOOM 

6.00 

DOOM 

15.00 



Dataset 

Data 

Missile 2 





AUE 101 

• 

o 





AUE 102 

■ 

□ 





AUE 105 

♦ 

o 





AUE 107 

▲ 

A 





Figure 6.- Continued. 


PUE 106 

CfiLSPflN CP PRD CO NTROL, 
B1M1C2T2 
BETQ 0.000 
PHI CM) 0.000 
PHI TflL 4S.C0C 
OOPOS 3.000 
TfEPOS 1.000 
DCH)1 0.000 

0002 -3.O00 

0003 0.D00 
DOOM -3.000 


ALE 109 

CfiLS°RN CRNW OONTROL. 
BIN1C2T2 
BETP 0.000 
PHIOO 0.000 
PHITPL 45.0C3 
0CP3S 3.0t0 
TfLPOS 1.000 
DOOl 0.000 
DOE2 0-.0O0 
0003 0.0C0 
DOOM 0.000 


PUE112 

CflLSPPN CRNflRD CONTROL, 
B1N1C2T2 
BETP 0.600 
PHICND 0.000 
PHITPL 45.000 
OCP03 3.000 
TPLPCS 1.000 
DODl 0.000 
0002 6.000 


PUE114 

CPLSPRN CfiNRRD CONTr.X, 
B1N1C2T2 
BETP 0.000 
PHIOD 0.000 
PHITRL 45.000 
OCPGS 3.000 
TflLFOS 1.000 
DOOl 0.000 
D002 15.000 

0003 0.000 

0004 15.000 


Dataset 
AUE 108 
AUE 109 
AUE 112 
AUE 114 


Data Missile 2 

• o 

■ o 

♦ O 

A A 

I 1 1 /' \ ' 

J \ 


A\ \ 






Canard Fins ^ 

Tail Fins X 

Moment Center - MS 32.5 
M_ * 0,9 


Angle of Attack, De<^ees 
(c) <f> * 0; interdigitated tail fins. 
Figure 8.- Continued. 
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Cft-SPfiN CP' .'PHD CONTROL, 


CflLSPAN CPK^D CONTROL, 
BIMC2T2 


CRLSPPN CflNAPD CONTROL, 
B1N1C2T2 


CflLSPff! CPNRRD CONVI 
B1N1C2T2 


BEIfl 

0 CO 

BETA 

0.C3 

BETA 

0.00 

BETA 

0.00 

PHIOO 45. CC 

PHIC-JD 

45.00 

PHIOO 

45.00 

PHIOO 

45.00 

PHITQL 45 ca 

PHITQL 

45. C0 

PHITRL 

45.00 

PHITflL 

45.00 

OC?OS 

3 00 

CNDF08 

3.C0 

OOPOS 

3.00 

ooros 

3.00 

TflLPOS 

1 00 

TflLPOS 

1.00 

TflLPOS 

1.00 

TflLPOS 

1.00 

DC? 01 

-3.00 

DOOt 

0.00 

DOC1 

6.00 

DOOl 

15.00 

DCJO? 

*3.00 

0002 

0.00 

0002 

6.00 

D0C2 

15.00 

DOD3 

-3.00 

0003 

0.00 

DOC3 

6.00 

DOOT 

15.00 

DOOM 

*3.00 

0004 

0.00 

DOOM 

6.00 

0004 

15.00 


Dr.taset 
AUE 115 
AUE 116 
AUE 119 
AUE 121 


Missile 2 


Canard Fins X 
Tail Fins ^ 

Moment Center • MS 32.5 


Angle of Attack, Degrees 
- 45°; interdigitated tail fins 
Figure 8.- Concluded. 






AUE006 

CEL SPAN CPNARD CONTROL, 
B1M1C2T3 
BETA 0.000 
PHICM) 0.000 
PHITAL 0.C00 
CNDPOS 3.000 
TALPOS i.eee 

0001 0.000 

0002 -3.000 

0003 0.000 
DOOM -3.000 


OUE088 

CAL£?AN CANARD QOhm, 
BIMC2T2 

BETA 0.000 
PHICM) 0.000 
PHITAL 0.000 
OCPOS 3.000 
TALPOS 1.000 
DOOl 0.000 
D002 0.000 

0003 0.000 
DOOM 0.000 


AUE081 

CALSPMN OANPH) CONTROL, 
81N1C2T2 
BETA 0.000 
MID© 0.000 
PHITAL 0.000 
CNJPOS 3.000 
TftPCS 1.000 
D>D1 0.000 
0002 6.000 
0003 0.000 
DOOM 6.000 


AUE09M 

CALSPR'. CANARD CONTROL, 
B1N1C272 
BETA 0.000 
PHICND 0.000 
PHITAL 0.000 
CM3P0S 3.000 
TALPOS 1.000 
DCfOl 0.000 

0002 15.000 

0003 0.000 
DOOM 15.000 


Dataaet Data 


Missile 2 



(a) $ ■ o; inline tail fins. 

Figure 9.- Comparison of pitching-moraent coefficient 
data and predictions from Program MISSILE 2 for 
Burt's Canard-Configuration tests; M* ■ 1.25. 


83 



PUE 101 

crlsppn cp'Ard control, 

B1N1C2T2 


BETA 

PHICND 

PHITflL 

OOPOS 

TflLPOS 

DOOl 

DCK02 

XK03 

doom 


e.00 

45.00 

0.00 

3.00 

1.00 
- 3.00 
- 3.00 
- 3.00 
- 3.00 


OE102 

CflLS^AN CRNhRD CONTROL. 
B1N1C2T2 


BETP 

PHIDO 

PHITQL 

OCPDS 

TftPOS 

DCND1 

D0O2 

DCN)3 

DCND4 


0.C3 

45.C3 

0.C3 

3.00 

1.00 
0.00 
0.00 
0.00 
0.00 


fluEtes 

CRLSPPN CflNCFD CONTROL, 
B1N1C2T2 


BETA 

PHICKD 

PHITRL 

OOPOS 

TflLPOS 

DOOl 

0002 

0003 

0004 


0.00 
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Figure 9.- Continued. 
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Figure 9.- Continued. 
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Figure 9.- Concluded. 
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r'igur* 15.- Flow chart for aubroutina CRKATO. 
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Figure 16.- Flow chart for subroutine FINDER 




Figure 17.- Flow chart for subroutine ABC. 
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